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ABSTRACT 

The ionizing background determines the ionization balance and the thermodynamics of the cosmic 
gas. It is therefore a fundamental ingredient to theoretical and empirical studies of both the IGM 
and galaxy formation. We present here a new calculation of its spectrum that satisfies the empirical 
constraints we recently obtained by combining state-of-the-art luminosity functions and intergalactic 
opacity measurements. 

In our preferred model, star-forming galaxies and quasars each contribute substantially to the HI 
ionizing field at z < 3, with galaxies rapidly overtaking quasars at higher redshifts as quasars be- 
come rarer. In addition to our fiducial model, we explore the physical dependences of the calculated 
background and clarify how recombination emission contributes to the ionization rates. We find that 
recombinations do not simply boost the ionization rates by the number of reemitted ionizing photons 
as many of these rapidly redshift below the ionization edges and have a distribution of energies. A sim- 
ple analytic model that captures the main effects seen in our numerical radiative transfer calculations 
is given. 

Finally, we discuss the effects of Hell reionization by quasars on both the spectrum of the ionizing 
background and on the thermal history of the IGM. In regions that have yet to be reionized, the spec- 
trum is expected to be almost completely suppressed immediately above 54.4 eV while a background of 
higher-energy (> 0.5 keV) photons permeates the entire universe owing to the frequency-dependence 
of the photoionization cross section. We provide an analytical model of the heat input during Hell 
reionization and its effects on the temperature-density relation. 

Subject headings: Cosmology: theory, diffuse radiation — galaxies: formation, evolution, high-rcdshift 
— quasars: absorption lines 



1. INTRODUCTION 

The cosmic baryons give the ultraviolet (UV) background a particularly important standing among radiation 
backgrounds. In fact, the ioniz ation potentials of b oth hydrogen and helium 3 , which together account for 99% of 
the baryonic mass density (e.g., |Burles et al.||200~lj ), correspond to electromagnetic wavelengths in the UV regime. 
The UV background therefore governs the ionization state of intergalactic gas and furthermore plays a key role in its 
ther mal evolution through photoheating. As such, it is an essential input to cosmological hydrodynamic simulations 
(e.g., Hernquist et al.|1996 Katz et al.|1996b Dave et al.|1999 Springel fe Hernquist|2003 ) as well as to observational 
studies of the intergalactic medium (IGM). 

The ionizing background can for example suppress the abundance of dwarf galaxies and the amount of cool 
gas in low-mass galaxies that do f orm both by modifying the cooling function through the ionization balance and b y 
heating the gas before it collapses ( Efstathiou|1992 Quinn et al.|1996 Thoul fc Weinberg|199Q _ Weinberg et al.|1997 l. 
It is also crucially important for any simulation ot the Ly a fores t, since the absence of aRjunn fc Peterson ( |1965 l 
trough in the spectra of quasars up to z ~ 6 (e.g., |Fan et al.||2002[ |2006b Becker et al.||2007p indicates that the IGM 
is highly ionized up to at l east that redshift. Since tFe optical depth of the Lyct forest is directly tied to the hydrogen 
photoionization rate fe.g., |Rauch et al.||1997] [McDonald fc Miralda-Escude|[200l| |Meiksin fc White||2003| |Tytler et al 
2004 Kirkman et al. 2005 Bolt on et al.||2005| |Faucher-Giguere et al.| 2008a b), it is important to know the latter 
accurately. In addition, the U V background determines the photoionization rate of helium, which is of particular 
relevance given the growing interest in studying Hell reionization, which may occur at redshifts z ~ 3 — 4 for which 
a wealth of observational data are already available and upcoming (<J7|. The full spectrum of the UV background is 
perhaps most important in the study of metal ions, such as SIV ana CIV, where relat ing the ionic abundances to 
elemental abundances or cosmic metal mass dens i ty requires ionizing corrections (e.g., Cowie et al. 1995 |Songaila| 
fc Cowie|[i"996l |Songaila| [200T| [Schaye et al.||2003] |Boksenberg et al.|[2003] |Aguirre et al.||2M4| |&mcoe et al-POMf 
Aguirre et al. 2007[ jFechner fc Richter 2008 ) . Finally, the spectrum of the UV background obviously depends on its 
sources and its st udy can therefore teach much about the sources responsible for keeping the IGM ionized, as well as 
reionizing it (e.g., |Miralda-Es cudc 20 03} |Bolton fc Hae hnelt 2007 Faucher-Giguere et al. 2008a|b[ ). 
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3 13.6, 24.6, and 54.6 eV for HI, Hel, and Hell, respectively. 
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Following early work (Miralda-Esc ude fc Ostriker 1990 |Shapiro et al. 1994 Giroux & Shapiro 1996[ ), Haardt 
& Madau ( 1996[ ) (see also |Fardal et al. 1998p pioneered calculations of the U V background spectrum in their study 
ol radiative transfer in a clu 



Haardt & Madau 2001 1 have since 



umpy universe. Their model and some variants (e.g. 
been extensively used in several hundreds of studies in the literature. Over a decade after their original calculation, 
the empirical constraints on the UV backgrounds and its sources have however improved dramatically. Larger and 



deeper surveys at all wavelengths have constrained the quasar luminosity function to both fainter magnitudes and 



higher redshifts (e.g., Boyle et al. 
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et al. 2004[ |Kichards et al.| 2005; 
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|Kichards et al. \ 
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2006; 


Hopkins et al. 2007 Fontanot et al.||2007 Siana et a 


|2008|). At the same time, our understanding of the popu 
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of high-redshiit star-forming galaxies has tremendousl y expanded thanks to the application of the Lyman break 



Yoshida et al. 



selection technique to ever more ambitious su rveys (e.g., |Steidel et al.||1999 Sawicki fc Thompson|2 006 
2006| Bouwens et al.|[2007 Reddy et al.| [2008). Detailed studies of the absorption properties of the IGM, particularly 
by HI and Hell, have also provided particularly valuable constraints on the UV background. These constraints are 



especially relevant for the UV background as the IGM is sensitive to the integral of the UV photons emitted by all 
sources, regardless of whether these are directly detected. Moreover, the IGM constraints probe the density of ionizing 
photons and thus circumvent the need to assume an escape fraction relating the luminosity of quasars and galaxies 
measured red ward of the Lyman limit to their net output of ionizing radiation. 



In a series of previous papers, we have measured the intergalactic Lya opacity (Faucher-Giguere et al. 2008d I 
and derived empirical constraints on the UV background and its so urces incorpo rating also information on the 
reionization of HI and Hell, as well as -/VHeii/-/VHi column density ratios ( |Faucher-Giguere et al.|[2008a|b[ ). Specifically, 
we found that the HI photoionization rate is remarkably constant over the redshift interval z = 2 — 4. Since the quasar 
luminosity function peaks strongly around z — 2, star-forming galaxies most likely dominate the ionizing background 
beyond z 3. The column density ratios however indicate that quasars likely do contribute a large fraction of the 
ionizing background at their peak. In this paper, we use these constraints as a basis for a new calculation of the full 
spectrum of the UV background. In addition to the improved empirical input, we reexamine many of the assumptions 
entering the original Haardt & Madau ( 1996 ) calculation. As we will show, we find that the original calculation likely 
overestimated the contribution of recombination emission to the ionizing background by a factor of a few. Of perhaps 
greatest interest, the original calculation completely neglecte d the effects of Hell reioni zation, so that simulators have 
usually resorted to artificial prescriptions to complement the Haardt & Madau ( 1996 1 spectrum. Here, we explicitly 
discuss the effects of Hell reionization on the UV background spectrum as well as on the thermal history of the IGM 
and provide a physical framework to implement them. 

We review the basic equations of cosmological radiative transfer and the column density distribution of HI ab- 
sorbers in Sj2] In Sj3] we study the ionization structure of individual absorbers and derive approximations to be used 
in the cosmological solution. £j4] is devoted to the calculation of the contribution recombinations to the cosmological 
cmissivity. Empirically calibrated calculations of the UV background spectrum, derived quantities, and their 
dependences on input parameters are presented in fj5] In Sj6] we investigate how the calculated spectra (including 
recombination emission) and the corresponding ionization rates depend on input parameters. The effects of Hell 
reionization are investigated in <|7j We finally compare our results with previous work in ^8] and conclude in EjOj 

A series of appendices supplement the main text with technical details. In Appendix [X] we describe our pho- 
toionization code. Appendix |B] contains technical aspects of our treatment of recombination emission, while Appendix 
[C] presents an analytic model of how this recombination emission boosts the photoionization rates. Appendix [D] 
analytically discusses spectral filtering in different regimes to aid in interpreting our results. Appendix [E] finally 
references atomic physics quantities used in our calculations. 

Throughout, we assume a cosmology with (Cl m , fib, ^a, h, as) — (0.28, 0.046, 0.72, 0.70, 0.82), as inferred 
from the Wilkinson Microw ave Anisotropy Probe (WMAP) five-year data in combination with baryon acoustic 
oscillations and supernovae (Komatsu et al. 2008). Unless otherwise stated, all error bars are la. Table [I] defines 
many symbols used here. 

2. COSMOLOGICAL RADIATIVE TRANSFER 

2.1. Radiative Transfer Equations 

In this work, we are first concerned with the specific intensity of the diffuse cosmological UV background averaged over 
both space and angl e, which we denote by J u . The basic equations of cosmological radiative transfer were particularly 
well summarized by Haardt & Madau (1996), on which we base our treatment below. The specific intensity satisfies 



the radiative transfer equation 



— - vH— 
dt dv 



J v = -3HJ U - cu v J v + — e„ 

47T 



(1) 
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TABLE 1 
Symbols used in this work 



Symbol 


Definition 


Tli 


number density of species i 


N, 


column density of species i 


Tv 


optical depth at frequency v 


f 


effective optical depth 


I V 


specific intensity along a ray 


J V 


angle-averaged specific intensity 


a/ 


absorption coefficient 


> 


emission coefficient 


T 


gas temperature 




photoionization cross section of species i 


r< 


photoionization rate of species i 


A,B 
a i 


case A or B recombination coefficient to species i 




recombination coefficient directly to level n = j of species i 


4>(u) 


line profile 


X 


mass fraction of hydrogen 


Y 


mass fraction of helium 


x u 


fraction of hydrogen in Hll 


yu, yui 


fractions of helium in Hcfl, Hclfl 



a The symbol a is also sometimes used as a recombination coefficient or 
a spectral index. The meaning should be transparent from the context. 

where H(t) is the Hubble parameter, c is the speed of light, a v is the proper absorption coefficient, and e„ is the proper 
emissivity. Integrating equation and expressing the result in terms of redshift gives 

J »o{ z o) = ^J dz^- ^"^3 e A z ) exp[-f(t/ , z , z)], (2) 

where v = vq(1 + z)/(l + zq), the proper line element dl/dz = c/[(l + z)H(z)], and the "effective optical depth" f 4 
quantifies the attenuation of photons of frequency isq at redshift zq that were emitted at redshift z by the relation 
e T = (e~ T ), where the average is over all lines of sights from z to z. For Poisson-distributed absorbers, each of column 
density iVni, 

f z r°° d 2 N 
f(v ,z ,z) = dz' dN H1 ——(l-e~^), (3) 
Jz Jo dN-sjdz 1 



where d 2 N / dNmdz' is the column density distribution versus redshift ( |Paresce et al~1l980| ). 



Note that these expressions neglect the clustering of sources and sinks of radiation in both the Poisson distribution 
assumption and in assuming that the spatial average (^(z) exp[— t(isq, z$, z)]) separates into (e„(z))(exp[— r(fo, zq, z)]) 
in the integrand of equation [2] The Poisson distribution assumption should be very good since the mean free path 
of ionizing photons of hundreds of comoving Mpc at most r edshifts of interest (see q7| far exceeds the correlation 
length < 5 comoving Mpc of the Lya forest absorbers (e.g., |McDonald et al. 2000 Faucher-Giguere et al. 2008c I 



While absorbers likely do cluster around sources, this effect can be viewed as being incorporated into the definition of 
the escape fraction. The above formalism will however break down in certain regimes where sources are rare and in 
particular during Hell reionization. We discuss these cases in ^j7] As in equation |2j we henceforth drop the explicit 
averaging brackets around the emissivity e u . 

The optical depth r„ shortward of the Lyman limit will be dominated by the photoelectric opacity of hydro- 
gen and helium, 

Ty = NHlO-HlM + ^HclCTHol(^) + NhcIICHcIi(^), (4) 

where the Ni and Oi are the column densities and photoionization cross sections of ion i. Only the distribution of 
Nhi is reasonably well determined over a large redshift interval. We will therefore make use of relations between iVjji 
and the column densities of helium established in 33] In our calculations, we will prescribe d 2 N/dNmdz as well as 



the specifi c emissivity of the ionizing sources, e s J c , based on our previous empirical studies (Faucher-Giguere et al. 
2008d|a|b| ) 



2.2. HI Column Density Distribution 

Following previous work and consistent with empirical constraints, we parameterize the column density distribution 
with power laws in TVhi and z: 

d 2 N _ / N Qilow N-f{l + z < z low 



dzdN m \N N-f(l + z)^ z>z low ' (5) 

4 This quantity is often denoted by r c ff. We use a different notation here to distinguish it from th e effective optical depth r c f{(z) = 
Tn(_F)(z) owing to Lya absorption measured from quasar spectra (e.g. , [Faucher-Giguere et al.|2008d| l. 
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The transition at zi ow accou nts for the flattening of the redshift evolution observed at z < 1.5 (e.g., Weymann et al. 
1998 Kim et al.||2001[ 2002ft th eoret ically under stood to arise from the d rop in intensity o f the ionizing background 
Snow redshifts f'l'heuns et al.||1998l |Dave et al.||l999| |Bianchi et al~1|2001| |Scott et al.||2002| ). We fix 7low = 0.2. 



For a steep column density distribution with ft < 2 
from systems of optical depth near unity (A/hi ~ 10 172 cm 
ft and 7 that are most appropriate in this neighborhood 
with C — 0.25 and 7 = 1.5 provides a good fit at leas t up to z 



to C by N = {ft- 



.min ' 



most of the contribution to r at the Lyman limit arises 
~ 2 ). We thus focus on the v alues of the power-law indices 
Stengler-Larrea et al.| ( |1995[ ) find that dN/dz = C(l + z) 7 
4. 1 for systems with iVm > 10 172 cm" 2 , whereas 



the column density power law is well-fitted by ft — 1.4 ( Misawa et al.|[2007 ). The constant N in equation [5] is related 



where Nui 



'HI 



10 



17.2 



cm . We use these values, with z\ a 



1.5, in fiducial 



calculations but explore varying these parameters in f|6] 



Before proceeding, we note that column density distribution is largely unconstrained above z — 4.1, the high- 
est redshift at which the Lyman limit system abundance has been measured. We therefore simply extrapolate from 
lower redshifts and caution that our calculation of the ionizing background may become inaccurate in this regime. In 
particular, we expect the extrapolation to become unreliable at z « 5.7, where the evolution of the effective optical 
depth m easured from the s pectra of z > 6 quasars diverges rapidly from the power law fitting the data below this 
redshift (|Fan et ah 2006b though see Becker et al. (2007) for an opposing point of view), perhaps owing to HI 



reionization. 

3. THE IONIZATION STRUCTURE OF INDIVIDUAL ABSORBERS 
3.1. Overview 

In this section, we study the photoionization equilibrium structure of individual cosmic absorbers composed of 
hydrogen and helium (with mass fraction 75% and 25%, respectively) as a function of the illuminating radiation 
background. This serves two purposes: close the set of equations of cosmological radiative transfer (fj2| and allow us 
to more realistically calculate the contribution of recombination lines to the ionizing background spectrum (^. To 
this end, we have developed a code that self-consistently solves the photoionization equilibrium balance, including the 
influence of recombination radiation. This code provides mor e accurate solutions than previous approximation s with 
semi-infinit e geometry and an escape probability formalism (Haardt & Madau 1996) or gray cross sections (i Fardal 
et al. 19981. To alleviate the text, the details of our photoionization calculations are provided in Appendix |A| We' 
assume our absorbers to be slabs of thickn ess equal to the Jeans scale of the gas, which is a function of t he assumed 
tempe rature T = 2x 10 4 K ( |Schaye|200l[ ). This temperature is consistent with the line-fitting analysis of [McDonald 
et al. (20011 and with the Lya forest power spectrum analysis of Zaldarriaga et al. (20011 for the z ~ 2 — 4 IGM at 
mean density. 



3.2. -/VhoII and Afjjei from AT] 



HI 



As only the column density distribution of HI is reasonably well constrained, the first application of our photoion- 
ization calculations is to obtain relations giving A'hcII and A?Hei m terms of A'hi- In Figure [TJ we show the numerical 



results for external spectra J„ 



10- 



erg s 



cm 



Hz~ 



sr 



>Mn)- 



ith a = 0.0, 0.5, 1.0, 1.5, 2.0, and 2.5 



from bottom up. To ensure that these test spectra are representative of the UV background, we suppress the power 
laws by a factor of 10 above the Hell ionization edge. 

Since these relations enter within three nested integrals (eqs. [TJ and Ejl, it is necessary to develop analytical 
approximations that are fast to evaluate. It would be impracticalto use self-consistent numerical photoionization 
calculations at each redshift and for each column density in the cosmological solution. Defining 77 = A^ieii/^HI, when 
both HI and Hell are optically thin and in the limit of nearly complete ionization we have 



?7thin = 



r He ii a& iX' 



(6) 



For fixed external background and increasing A^hi, an absorber first becomes optically thick in Hell, at which point 
77 increases rapidly with Nm- The absorber then becomes optically thick in HI as well and, owing to the greater 
abundance of hydrogen, ATjji finally rapidly overtakes AThcII- This leads to the plateau, increase, and then decrease 
of rj with respect to A'hi seen in the numerical calculatio ns. Similar behavior is found in three-dimensional radiative 
transfer simulations of the IGM ( Maselli & Ferrara 2005 ) . 



Fardal et al. (19981 give a fitting formula derived under the assumptions of negligible AThcI and nni/nn <C 1: 

Y 



16X1 + At hi 



Im — r Hcii 



THell 



1 + BthcII 



Hell, 



(7) 



where I\ = r? xt /rt e a^, thi = ohiAhi and t Hell = CHeiiA^Hcii, and we have generalized their result to allow for 
arbitrary coefficients A and B. These fitting coefficients depends on, in particular, the relation between A^hi and 
n e and we will not be using the same model as these authors. Although our numerical calculations do not a priori 
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assume a relation between iVui and n e , we must assume one in order to make use of the analytic approximation in 
equation[7j The approximation curves in Figure [ljcan be reproduced by taking n e = 1.4 x 1CP 3 cm -3 (7Vhi/10 17 ' 2 
cm~ 2 ) 2 / 3 (THi/10~ 12 s -1 ) 2 / 3 . This relation is approximately derived under the assumptions of Jeans length thickness 
of the absorbers and optically thin photoionization equilibrium at T = 2 x 10 4 K. Figure [l] shows that A — 0.15 and 
B = 0.2 give a good fit to our numerical results for a wide range of external illuminating spectra. The fitting formula 
has the exact optically thin limit; the asymptotic divergence from the numerical results as -ZVhi — > oo is unimportant 
as most of the Hell opacity arises in systems with thi < 1. Although we have as sumed specific (but varying) spectral 



shapes in determining the fitting parameters A and B, the Fardal et al. (1998) derivation of the functional form in 
equation [7] illustrates how the relation between A^jeii and iVni depends principally on the photoionization rates Thi 
and ThcH- The relation should therefore hold well in general. 

Obtaining a physically- motivated analytic approximation to £ = Nn e i/Nm is more difficult since yucH is n °t 
readily known (for 77, we know that xhii ~ ym ~ 1 in the almost-complctcly ionized case, so that the ionized fractions 
do not appear explicitly in eqs [6] and [7]) . Because the ionization potential of Hel is relatively close to that of HI, their 
ionization states are similar and since helium is less abundant by a factor of 12 by number, Hel should contribute 
relatively little to the ionizing opacity. This intuition is supported by the right panel of Figure [T] which shows that £ 
versus Nm is 1 for illuminating spectra considered. After hardening by IGM filtering, both star-forming galaxies 
and quasars are expected to produce roughly flat spectra between the HI and Hel ionization edges (fj5j), so that a 
representative case is a « 0, yielding £ < 10 -3 . We therefore approximate £ = in our cosmological calculations and 
verify using toy cases of constant £ that this is justified in £j6] 




Fig. 1. — Ratio r\ = -/VhoIi/JVhi (left) and £ = Nn c i/Nm (right) as a function of HI column density. The solid curves show full numerical 
photoionization calculations and the dashed ones show analytical approximations based on equationpjfor r). From the bottom up, power-law 
external spectra J„ = 10 — 21 erg s _1 cm" " 2 Hz" 1 ST- 1 (i//v m )~ a with a = 0.0, 0.5, 1.0, 1.5, 2.0, and 2.5, suppressed by a factor of 10 
above the Hell ionization edge, are assumed. 



4. RECOMBINATION EMISSION 
4.1. Cosmological Emissivity 

The cosmic absorbers not only act as sinks but also as sources of ionizing radiation as a certain fraction of 



& Madau 


1996 


Fardal et al. 


1998 



This recombination emission must be taken into account because it may boost 
the photoionization rates and also since line reemission can imprint significant narrow features in the ionizing 
background spectrum. Our approach to include this recombination contribution is based on the self-consistent 
numerical calculations of recombination emission from indi vidual absorbers using th e code outlined in the previous 



section and detailed in Appendix [A] This again differs from |Haardt fc Madau ( 1996| , who used an analytical escape 
probability formalism and made the assumption of a constant sourc e function wi t hin th e absorbers (which breaks 
down in the very optically thick systems), and from the treatment of Fardal et al. (1998) and thus provides a check 
of these results. 



For each recombination process of interest, we calculate the emergent specific intensity IJ ec (A r Hi) owing to the 
process given the external illuminating spectrum numerically using our photoionization code. The cosmological 



G 



Faucher-Giguere et al. 



recombination emissivity for this process is then an average over the column density distribution: 

^ c = 4 ^r^ Hi ^ /r(7vm) - (8) 

Since in general /™ c (./Vhi) depends on the spectrum of the ionizing background, which is not known a priori and 
evolves at each step in the redshift integration, it is again necessary to obtain an analytical approximation for this 
function that scales appropriately with the external background, as it is not practical to perform self-consistent 
numerical calculations in the cosmological solution. We develop these analytical approximations in Appendix [B] for 
each recombination process of interest. 




Fig. 2. — Ratio of the HI photoionization rate outside an individual absorber that is contributed by different recombination pro- 
cesses (HI LyC in blue, Hell BalC in red, Hell Lyct in magenta, Hell LyC in green) to the external photoionization rate. The 
solid curves show full numerical integ ratio ns over photo-ionized slabs (eq. |B2[ | and the dashed one s sho w analytical approximations 
based on the optically thin limit (eq. |B3| and saturation in the optically thick regime described in j ]B.2| External power-law spectra 
= 10 — 21 erg s -1 cm -2 Hz -1 sr _1 (y/i/pjl) - a with a = 1.0, 1.5, 2.0, and 2.5, suppressed by a factor of 10 above the Hell ionization 
edge, are assumed in the different panels. 



4.2. Recombination Processes 

For hydrogen, the only ionizing process is direct recombination to the ground state, which produces a 1 Ryd HI 
LyC photon. For helium, both recombinations to Hell and to Hel can in principle p roduce ionizing photons. In 



cosmological conditions, Hel plays a negligible role (j |3.2| and [6] Haardt & Madau 1996 1 and we will ignore it in our 



reemission calculations. Three permitted Hell recombination channels lead to the reemission of ionizing photons: 
Hell LyC recombinations directly to the ground state, indirect recombinations leading to Hell Lya emission, and 
recombinations to the n = 2 excited level resulting in Balmer continuum (BalC) photons. These respectively give 
photons of energy 4, 3, and 1 Ryd. Higher Hell Lyman-series photons could also produce ionizing photons, but we 
assume case B conditions in which these are ultimately degraded into lower-energy photons, the only ones of which 
that can ionize hydrogen being Hell Lya. We do not include forbidden two-photon recombination processes as these 
are energetically subdominant and do not result in distinctive emission features. 



In calculating the contribution of reemission to the photoionization rates, it is important to model the finite 
width of the recombination lines. If HI LyC reemission is incorrectly modeled as a 6— function, the reemission photons 
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Fig. 3. — Frequency integral of HI Lya reemission as function of HI column density . Th e solid curves show full numerical photoionization 
calculations and the dashed ones show analytical approximations based on equation |B8| From the bottom up, power-law external spectra 
= 10 -21 erg s — 1 cm -2 Hz -1 sr — 1 (v/i'hi)~ a with a = 0.0, 0.5, 1.0, 1.5, 2.0, and 2.5, suppressed by a factor of 10 above the Hell 
ionization edge, are assumed. 

are immediately redshifted below the HI ionization edge and are lost as contributors to the ionizing background. For 
continuum recombinations, the line profile is well approximated by 

, , \ (^/^ 10 c) _1 exp(-/i^/fcr) B(v - V Iec ) 

Precis) = r 7T7f\ ) (9) 

r(0, hv rcc /kT) v rec 

where T is the temperature of the gas and 9(Av) is the Heaviside function which is 1 for Av > and otherwise 
(Appendix [E| . Electrons in gas at higher temperature tend to have large kinetic energy and so give rise to higher- 
energy recombination photons that take longer to redshift below the ionization edge. In Appendix [E] we show that 
broadening owing to thermal and peculiar motion is negligible relative the width of the profile in equation [9] For 
Lya emission, either by HI or Hell, a 5— function profile </>hi/HcII Lya^) = 5{v — t'Hi/Heii Lya) is however appropriate 
because of the narrow intrinsic line width (much smaller than the mean free path by which photons are redshifted 
before being reabsorbed) and its distance from the ionization edg es. While resonant scattering radiative transfer 



effects can broaden Lya emission lines b y - 10- 1000 km s^ 1 (e.g., |Neufeld||l990| |Zheng fc Miralda-Escude|2002~ 



Di- 



jkstra et al.|2006| Verhamme et al. 2006 1, this width is negligible in comparison to the cosmological redshift broadening 



The analytical approximations for the recombination emission from individual absorbers are compared to the 
full numerical calculations in Figure [2] for the ionizing processes. In all cases, r^/r^' (where we define 
^Hi(-^Hi) = 47r dv /{hv)I r v cc {NY{i)<Jm{v)) is maximum for HI LyC reemission, as expected since hydrogen 
recombinations are more frequent owing to its greater abundance and these recombination photons have the largest 
photoionization cross section, and is equal to about 10%. The helium recombination processes all contribute at the 
10~ 3 level or less. Note, however, that Hell LyC reemission will contribute more significantly to the Hell ionizing 
background and that processes which contribute negligibly to the photoionization rates can still imprint important 
narrow features in the background spectrum that can be important for metal line studies. The agreement between the 
numerical calculations and analytical approximations is generally good and the approximations scale well for different 
spectral indices. Discrepancies of a factor of a few exist over some column density intervals, particularly for the Hell 
BalC and Hell LyC processes. These processes are complex in their details that depend on the non-monotonic relative 
ionization of hydrogen and helium (Fig. [lj but their contributions are nevertheless reasonably well captured and 
fortunately subdominant to the photoionization rates. In contrast, the dominant contribution of reemission to the 
hydrogen ionizing background, HI LyC emission, involves only hydrogen and is accurately and robustly approximated. 

Figure [3] compares the analytical approximation for HI Lya approximation to the full numerical solutions. In 
this case, both the optically thin and optically thick limits are accurately captured, resulting in an excellent 
approximation at all column densities that scales correctly with the external illuminating spectrum. 



5. EMPIRICALLY CALIBRATED SPECTRA 
5.1. Quasar and Stellar Emissivities 



Having established efficient approximations for the radiative transfer within individual absorbers (^J3]and|4|, we 
proceed to include these in the solution of the cosmological radiative transfer pro blem (jj2j). Our prescrip tions for the 
sources of ionizing radiation are based on the empirical constraints obtained in Faucher-Giguere et al.| ( |2008d|aTr7 



Note, however, that these prescriptions can easily be modified to accommodate further constraints: our numerical code 
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can compute the ionizing background for arbitrary input emissivities. We explore variations about these fiducial pa- 
rameters in Ej6j Here, we consider two dominant known sources of ionizing radiation: quasars and star-forming galaxies. 



For the quasar emissivity, we use the quasar luminosity function of Hopkins et al. ( 2007 ) based on a large set 
of observed quasar luminosity functions in the infrared, optical, soft and hard X-rays, as well as emission line 
measurements. Denoting by e# the emissivity at 4400 A and assuming Lb = vL v \ 



'4400 A 



QSO,com 



dL f 



d(f> Lb 
cLLb v\ 



(10) 



4400 A 



where d(f>/dLB is the £?-band luminosity function in comoving units. The emissivity shortwa rd of 4400 A is calcu lated 
assuming that quasars have a spectral index a = 0.3 at 2500-4400 A, 0.8 at 1050-2500 A dMadau et al.||l999l), and 



chqso shortward of 1050 A. In order to match the total HI photoionization rate measured from the Lya forest and 
to account for uncertainties in converting from the e miss ivity at 4400 A to the photoionizati on rate, we allow this 
emissivity to be normalized by a constant factor (see S5.2I. In our fiducial model, aQso = 1-6 ( Telfer et al.||2002 l but 
note that other studies have found both softer and harder spectra, with a significant variance about the mean (e.g., 
Zheng et al.|l997| |Scott et al.||2004[ ). 



For the stellar emissivity, we assume that the emissivity is proportional to the star formation rate density, 

,t.com Tv- -com 

with the observationally-calibrated proportionality constant accounti ng for the efficiency of c onver sion of mass into 
ionizing photons. 



(11) 



We use the theoretical star formation history of Hcrnquist & Springel 

2003J) and simple analytical arguments 



a combination of hydr o dynam ical simulations (Springel & Hernquist 
Faucher-Giguere et al. (2008b), we found that this model provides a 
the Lya forest over 



2003]) developed from 
In 



setter fit at high redshifts to the opacity of 
1.2, is easier to reconcile with hydrogen reionization completing by z = 6, and is in 



z = 2" 

better agreement with the rate of long gamma-ray bursts observed by Swift than many of the existing measurements 
based on galaxy surveys, among which there is still a wide dispersion. We use equation 45 of |Hernquist fc Springel 
(20031 to scale their fiducial model t o the WMAP5 cosm ology assumed in this work. We assume, motivated by the 
theoretical starburst calculations of Kewley et al. (20011, that star- forming galaxies have a spectral index a* = 1 
bet ween 1 and 4 Ryd. This mod el is applicable lor the stellar populations calculated with the PEGASE code using 
the Clegg & Middlemass (19871 atmosphere models for Wolf-Rayet stars. While different theoretical assumptions 
lead to significant variance in the 1—4 Ryd spectrum, t his model prov i des th e best observational match to the hard 
starburst spectra inferred by optical line diagnostics by Kewley et al. (20011. We assume that they effectively emit 
no harder photons, the theoretical calculations showing a break of several orders of magnitude at the Hell ionization 
edge. In the UV spectrum redward of 1 Ryd, we take a* = 0, consistent with the LBGs observed by Shaplcy et al. 
(2003). Finally, we assume that the stellar emissivity has a discontinuity of a factor of 4 at the Lyman limit. While 
this factor is neither well constrained empirically or observationally, it only affects our predicted spectra (normalized 
to the measured ionizing background) at energies less than 1 Ryd, which we do not attempt to accurately model in 
this work. 

The emissivities are converted to proper units before being inserted in the solution to the cosmological radia- 
tive transfer solution in equation [2] and the total emissivity is then e„(z) = e^} so (z) + ej(z) + e£, cc (z). 

5.2. Results 

In Figure [4j we show the calculated cosmological UV background spectra at z — — 5 for the fiducial model above, 
with and without the recombination processes included. Since only quasars are assumed to contribute photons above 
the 4 Ryd Hell ionization edge, only them contribute to the Hell recombination lines and photoionization rate. 

In Figure [5] we show the integrated photoionization rates of HI and Hell, as well as the fractional contribu- 
tion of recombination lines with respect to the total background including both stars and quasars. The quasar 
contribution to the HI ionizing background increases toward z ~ 2 as the peak of the quasar luminosity function 
is approached; the z > 3 photoionization rate is dominated by stellar emission. The fractional recombination 
contribution to the HI photoionization rate ranges from 5% to 17% over the interval z = — 6, significantly smaller 
than the aHi,n=i(T)/ a Hi(^ 1 )lT=20,ooo k ~ 48% fraction of HI recombinations that are directly to the ground state. 
The relatively small contribution of recombinations to the ionizing background owes to a combination of the saturation 
of reemission in optically thick systems (Fig. j2J), leakage of the reemitted photons at the ionizing edge, and the 
frequency dependence of the photoionization cross section 6.2 and Appendix O). 



The fractional contribution of Hell recombinations to the Hell photoionization rate is also relatively small for 
the same reason, but is more difficult to calculate accurately at redshifts z > 4 in our model. In order to obtain 
an accurate value, in addition to the Hell LyC reemission line to be well resolved on the computational frequency 
grid, the mean free path of Hell ionizing photons must also be well resolved by the redshift grid. In our calcu- 
lation, quasars produce a negligible and rapidly dropping Hell photoionization rate at z > 4 while star-forming 
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Fig. 4. — Spectra of the UV background obtained by solving the cosmological radiative transfer equation (cq. rfj at different rcdshifts 
for our empirically-calibrated fiducial model with star-forming galaxies and quasars described in jj5] The thin curves ignore recombination 
emission by intergalactic absorbers and the thicker curves include this contribution. Star-forming galaxies dominate the HI photoionization 
rate at z > 3, with the quasar contribution becoming more important as the z ~ 2 peak of the quasar luminosity function is approached. 
Only quasars are assumed to produce Hell ionizing photons and only they contribute to Hell recombination processes. The integrated 
photoionization rates are given in Figure [5] 

galaxies maintain a roughly constant HI photoionization rate. In these conditions, the ratio rj = Nn e ii/Nffl tends 
to infinity and the Hell mean free path to zero, making it exceedingly difficult to resolve it. Fortunately, the 
total Hell photoionization rate in this regime is so small that its fractional enhancement from recombinations 
is of little practical importance. Moreover, in this regime Hell rcionization may well be still underway and the 
Hell ionizing background consequently modified, as elaborated on in SjT] In Figure [5] we indicate the portion 
of poor convergence by a dashed curve segment; the turnover of rggjj ICC /rj!f° I j cc around z ~ 5 is likely an art i- 
fact and we in fact expect it to continue to increase slightly toward higher redshifts owing to the reduced leakage (^|6.2[). 



est at z 



4.2, 



The total HI photoionization rate matches the value Thi = (0.5 ± 0.1) x 10~ 12 s _1 derived from the Lya for- 
subject to the constraint that quasars must contribute a large fractio n ne ar their peak 
This was done by normalizing the nominal quasar contribution (j |5.1| ) by a factor of 
:ontribution so as to provide the rest of the ionizing photons. The renormalization 
of the quasar contribution can be justified by uncertainties in the mean free path of HI ionizing photons (a direct 
product of the pre scribed HI column density distr ibution) , in their escape fraction, and in the quasar spectral template 



(Faucher-Giguere et al. 


2008a b) 


0.36 and normalizing t 


le stellar 



(see discussion in Faucher-Giguere ct al 
luminosity to the photoionization rate 



2008b) 



These uncertain factors enter in the conversion from the quasar 
"Since we wish to reproduce the more robustly constrained photoionization 
rate measured from the Lya forest, we adjust the normalization accordingly. 



Although our calculations are normalized to match the hydrogen photoionization rate measured from the Lya forest, 
it is important to emphasize that this measurement and hence the normalization of the spectr a calculated here ar e 



somewhat uncertain. The measurement was obtained using the flux decrement method (e.g., [Rauch et al. 

of tfe 



1997). 

— Tl — 

1 wo 



in which we solve for the value of Thi needed to produce the measured mean transmission of the Lya forest 
important sources of systematic uncertainty are the assumed IGM temperature (since the flux decrement constrains 
only the combination Thi /a^ T (T)) and the gas density distribution (whose details depend on the cosmological 
parameters and thermal history). Another potential worry is that the measured Lya forest mean transmission may 
be increasingly biased high toward high redshifts (inducing a redshift-dependent error) as the continuum level is 
increasingly absor bed and difficult to estimate d irectly. We have however quantified and corrected for this effect in 
our measurement (Faucher-Giguere et al. 2008d I and so it should not affect our results. In the end, we expect the 
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Fig. 5. — Photoionization rates for the UV background shown in Figure]?] Top left: HI photoionization rates including the contribution 
from recom binations. The total (quasars+ stars) photoionization rate is compared to the value inferred from the Lya forest forest flux 
decrement | |Faucher-Giguere et a l. 2008a b). Top right: ratio of the total HI photoionization rate to the value value obtained neglecting 
recombination emission] the bottom panels arc analogous but for Hell ionizing radiation. In this case, quasars are t he only contributors. 
The dashed portion of the rgg t I j roc /r^ [ ° r I j c curve at z > 4 indicates a regime of poor numerical convergence (see [5.2 I and the turnover at 
z ~ 5 is likely an artifact. 

measured Thi to be accurate within a factor ~ 2, with the possible errors mostly systematic and weakly dependent on 
redshift. For a more exhaustive discussion of the uncertainties of the measured Thi, see §3 of Faucher-Giguere et al. 



(2008b) 



Finally, we must also note that for precise work with hydrodynamical simulations, the simulated Lya forest 
mean transmission should always be compared with the measured value and the photoionization rates renormalized if 
necessary. In fact, even if the correct ionizing background (with the correct normalization) is prescribed, the simulated 
mean transmission may be slightly off if, for example, the temperature of the IGM is incorrect reproduced. This is 
particularly likely to occur if the effects of HI and Hell reionization (see Q are not explicitly modeled. 

6. DEPENDENCES ON INPUT PARAMETERS 

Even after fixing the stellar and quasar emissivities for our fiducial model, the spectrum calculations depend on 
a number of parameters. It is useful to investigate how the calculated spectrum and its integrals depend on these 
as their values are only known to limited precision. This also provides a physical understanding of the shape of the 
calculated spectra. We begin by considering the dependences of the overall ionizing background spectrum in §6.1| and 
focus on the recombination contribution in |6.2| 
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Fig. 6. — Effects of parameters on the calculated spectrum. Top left: varying the constant ratio £ = Vhci/Vhi = 
0, 0.001, 0.005, 0.01, 0.05, and 0.1. Top right: varying the HI column density distribution power-law index j3 = 
1.2, 1.3, 1.4, 1.5, 1.6 and 1.7. Bottom left: varying the stellar spectral index a* = 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0. Bottom 
right: varying the quasar spectral inde x a p ro = 0.5. 1.0, 1.5, and 2.0. In each panel, the curves correspond to these values from the top to 
the bottom. For realistic values of ( (j |3.2| and Fig. [lj for a background spectrum arising from stars and quasars, the effect of Hel is small. 
The column density distribution power law /3 determines the spectral hardening just shortward of the ionization edges. The stellar and 
quasar spectral indices determine the spectral slope of the background. For fixed emissivity at the Lyman limit, the stellar spectral index 
has only a modest effect on the amplitude of the spectrum because it is truncated at 4 Ryd. The quasar spectral index, assumed to extend 
to infinity, has a more drastic overall impact toward high energies. Recombination emission has been omitted for clarity of presentation. 

6.1. Overall Spectrum 

In Figure [6] we show how the spectrum changes when the constant ratio £ = -/Vh i/-/Vhi, the HI column density 
distribution power-law index [3, and the stellar and quasar spectral indices a* and chqso ar( 3 individually varied. In 
each case, all other parameters are fixed to the fiducial model of the previous section. Even for a constant ratio 
( = 0.1, a factor more than one hundred times that expected in our fiducial calculation (j |3.2| |, Hel absorbs only a 
very small fraction of the spectrum shortward of its ionization edge. It is therefore a good approximation to neglect 
it in our cosmological calculations. The HI column density distribution power-law index (3 determines the spectral 
hardening just above the ionization edges following a — > a — 3(j3 — 1) (Appendix |D| as well as the depth of the 
absorption edges. Note that the depth of the Hell absorption edge is more sensitive to (3\ th is arises because the 
column density distribution is normalized to the abundance of HI Lyman limit systems (j |2.2[ ) so that it is fixed in 
these calculations while the abundance of the Hell Lyman limit systems varies. The stellar and quasar spectral 
indices simply determine the spectral slopes of the background prior to hardening. For fixed emissivity at the Lyman 
limit, the stellar spectral index has only a modest effect on the amplitude of the spectrum because it is truncated 
at 4 Ryd. The quasar spectral index, assumed to extend to infinity, has a more drastic overall impact toward high 
energies: as v — > oo and spectral hardening becomes negligible, different spectral indices result in a J v ratio of 
(zv/^Hi) aQSO ' 1-QQSO ' 2 • At 10 keV, for example, this ratio is 735 for aQso.i = 1-5 and aQso,2 = 0.5; redshifted from 
z = 2 to z = 0, this falls in the bandpass of x-ray observatories such as C handra and XMM-Ne wton. As most (> 80%) 
of the soft x-ray background has already been resolved into AGNs (e.g., |Hasinger et al.||2005[ ), the x-ray background 
is a powerful probe of the high-energy quasar spectral energy distribution, although a proper analysis requires the 




Fig. 7. — Effects of the redshift evolution of the HI column density distribution on the hydrogen photoionization rate. Left: varying the 
redshift at which the column density distribution transitions to relatively flat evolution zi ow = 0, 0.5, 1.0, 1.5, 2.0 (bottom up). Right: 
varying the power-law index of redshift evolution 7 = 1.0, 1.25, 1.5, 1.75, and 2.0 (top-bottom) of the column density distribution at 
z > 2 low Varying z\ avl at z < 2 has little impact on Thi since at low redshift the HI ionizing mean free path is sufficiently large that the 
spectral intensity is limited by the cosmological horizon. The high-redshift Thi declines more rapidly with more a rapid increase in the 
abundance of absorbers with redshift, or large 7, translating into a more rapidly diminishing mean free path. 



inclusion of obscured quasars, which we do not explicitly consider in this work (e.g., |Gilli etld~1[2007t . 



In Figure [7j we explore how the hydrogen photoionization rate is affected by the redshift evolution of the col- 
umn density distribut ion. In the left panel, we vary the redshift at which the redshift evolution of the column density 
distribution flattens ( p.2| from z\ ow = to z\ OVJ = 2. Interestingly, this has a minimal impact on the redshift evolution 
of the photoionization rate even if it does significantly change the mean free path of HI ionizing photons at these 
redshifts. This is easily understood as a consequence o f the fact that the universe effectively becomes transparent at 
a "breakthrough" redshift Zbt ~ 2 ( |Madau et al.|[l999 l, below which the mean free path becomes so large the local 
ionizing background is not limited by the latter but by the cosmological horizon. As shown in the right panel, the 
high-redshift Thi declines more rapidly with more a rapid increase in the abundance of absorbers with redshift, or 
large 7, translating into a more rapidly diminishing mean free pa th. At present, although more than a decade old, the 



1994 1996 Stengler-Larrea 



best constr aints on the abundance of the Lyman limit systems (|Storrie-Lombardi et al 
et al.|[l995 l are relatively loose and mostly nonexistent beyond z 
toward higher redshifts, it is possible that these will give more credence to one of the alternative values of 7 plotted 
here. 



4. As future measurements refine these and push 



6.2. Recombination Contribution 

The contribution of recombinations to the photoionization rates, r wlth rec /r no rcc ; is a subtle question as it depends 
on several factors. It not only depends on the number of rcemitted ionizing photons integrated over the distribution 
of absorbers (Q but also crucially on the energy at which these photons are reemitted as well as on their redshifted 
energy at the point of evaluation of the photoionization rate. 

The LyC recombination line processes, most important for the boosting the ionization rates, reemit ionizing 
photons just above the ionization edges of HI or Hell. Since the ionizing background at a given point is sourced 
along its past light cone, its photons have generally redshifted slightly from their emission energy. As a result, many 
recombination photons with initial energy just above their corresponding ionization edges quickly redshift below these 
edges and are lost as contributors to the ionization rates. The fraction of ionizing recombination photons lost in 
this way depends on two factors: 1) the recombination line profile which determines how far above the ionization 
edge an ionizing photon is reemitted and 2) the mean free path of ionizing photons which determines how long the 
photons have to redshift before they are absorbed. In the limit of a mean free path of zero length, the recombination 
photons cannot redshift before they are reabsorbed and no photons are lost; as the mean free path increases, more 
photons leak out of the ionizing range. Similarly, a narrow line profile concentrates the recombination photons just 
above the ionization edges, leading to a high probability of leakage, while a wider one allows them to remain longer 
in the ionizing range. Since the mean free path of ionizing photons is determined by the column density distribution 
and the recombination line width is determined by the temperature of the absorbers, these are important parameters 
for the recombination contribution. Finally, as recombination photons are reemitted at energies above the ionization 
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Fig. 8. — Dependences of the recombination contribution on HI photoionization rate. Shown are ratios of the total HI photoionization 
rate, including recombination emission, to the same calculation ignoring recombination emission. In all cases, the sources of the ionizing 
background are fixed to the fiducial model of j}5] but we vary the parameters of the column density distribution and the temperature of the 
absorbers. Top left: C = 0.125, 0.25, 0.5, 0.75, and 1 from bottom up. Top right: = 1.2, 1. 3. 1.4. 1.5. 1 
Bottom left: T = 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 X 10 4 K from bottom up. Bottom right: |Haardt fc M; 
distribution with different redshift evolutions in the optically thin and optically thick regimes. 

edges and subsequently redshift, the frequency dependence o; L (y) oc v~ 3 of the photoionization cross section changes 
the weight they receive in the photoionization rates. 

All these effects are self-consistently treated when solving the radiative transfer equation [T Figure [8] shows 
how much the HI photoionization rate is increased by recombination emission as a function of the normalization C 
of the column density distribution, its power-law slope (3 (see eq. pL and the temperature T of the absorbers. We 
also show the case of the column density distribution assumed by [Haardt fc Madau ( 1996 ) , in which optically thin 
and optically thick absorbers have different redshift evolutions, leading to a redshift-dependent effective slope of the 
distribution (steeper at high redshifts). 

Because the mean free path decreases with increasing abundance of Lyman limit systems (the normalization 
C), fewer recombination photons leak out of the ionizing range and so the ratio of the photoionization rate with and 
without recombinations, Tgj th rec /rff T rec , increases. The ratio also increases with the steepness of the column density 
distribution since more recombinations occur in optically thin systems, from which practically all the recombination 
photons escape into into the IGM, as opposed to in optically thick systems that trap a large fraction. Higher gas 
temperatures result in wider recombination lines so that fewer photons are lost owing to redshifting as well as a 
higher fraction of recombinations directly to the ground state (aHi.n=i(T)/aHi a=i(-^)j Appendix |cj) . A competing 
effect is that the recombination photons of high-temperature absorbers tend to have higher energies and receive less 
weight in the photoionization rate. The net effect is however relatively weak on rgj th roc /rj!g rcc for the relevant 
temperature T ~ 2.0 x 10 4 K. Note that for any given set of parameters the ratio tends to increase toward higher 
redshifts since the mean free path is lower at the higher cosmological densities. This behavior is however not seen 
below the breakthrough redshift Zbt ^ 2, where the photoionization rate is limited by the cosmological horizon rather 
than by the mean free path. 
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In Appendix [C] wc develop a quantitative analytic model that captures and clarifies these effects and agrees 
well with the full numerical calculations presented here. 

7. REIONIZATION EVENTS 

The previous calculations have implicitly assumed that the universe is reionized in both HI and Hell. This 
assumption is most evident in the case of HI, for which we have used a column density distribution measured from the 
z < 4 Lya forest. The assumption creeps in for Hell rcionization during which there are large Hell patches the inside 
which the Hell photoionization rate is very low in comparison to within ionized bubbles. In each region, the mapping 
between iVjji and A^Heii depends on the local spectrum and will in general be very inhomogeneous. The IGM opacity 
to Hell ionizing photons will therefore be poorly approximated by using a globally-averaged spectrum to map from 
Nm to AThcii- In this work, we do not attempt to model HI reionization (an d the likely simultaneous reionization 
of Hel), which occurs at the limit of the present observational reach at z > 6 (e.g., Fan et al. 2006a Dunkley et al. 2008 1. 



The reionization of Hell was however likely delayed until the rise of the quasar luminosity function, at red- 
shifts that are immediately accessible to observations and some understanding of its effects can be obtained by 
studying the ratio r\ = -ZVhoIi/-^hi- ^ n f act > while stellar spectra are theoretically expected to have a strong break at 
the Hell ionization edge and therefore hav e litt le impact on the Hell ionization state, quasars have power-law far-UV 
spectra that extend well into the x-rays ( §5T| . Theoretical c alculations based on the quasar luminosity function in 
fact indicate that quasars can reion izc Hell by z <~ 3 — 4 (e.g., Sok asian et al.||2002| |Wyithe fc Loeb||2003| |Furlanetto" 



& Oh 2008 F aucher-Giguere et al.|[2008b[ ) . A number of lines of evidence, based HI and Hell Lya forests as well as 
on the evolution of metal line ratios, also suggest that the IGM is undergoing changes that could be associat ed with 



Hell reionization at these redshifts (for a review of these lines of evidence, see Faucher-Giguere et al. 2008d I . While 
alternative candidate s ources of Hell reionizatio n exist - such as possible H ell ionizing emission from galaxies (e.g., 
|Furlanetto fc O h 2008), h igh-re dshift x-rays (e.g., Oh|2001b Ricotti fc Ostriker|2004 ), or thermal emission from shock 
heated gas QMiniati et al.||2004[ ) - quasars are the best established and most likely. Large fluctuations observed in the 
Hell ionizing background toward z = 3 , which can be explained by the small number density of bright objects, also 



lend support to the quasar hypothesis ( Zheng et al. 2004 Shull et aL"||2004 Bolton et al.||2006 1 and we will therefore 
concentrate on this scenario. 



7.1. The Ionizing Background During Hell Reionization 



Recently, McQuinn et al.| ( 2008 1 performed detailed radiative transfer simulations of Hell reioni zation in large 
boxes up to 430 comoving Mpc on a side (for previous simpler tre atments, analytic and numerical, seelSokasian et aT 



Gleser et al. 2005 Paschos et al. 2007 Furlanetto & Oh 



2008 



Bolton et al. 2008 



These 



[20021 |Bolton et al.||2004| 

simulations used realistic models for the quasar sources based on the luminosity function of Hopkins et al. (20071 



and with physically and empiri cally motivated prescriptions for the triggering of quasars in massive halos as well as 



of quasar light curves (see, e.g., Hopkins et al. 2005a|b 2006 2008). A striking result of this work is the remarkable 



complexity of Hell reionization, in particular of the Hell ionizing radiation field, likely rendering the detailed resulting 
structure beyond analytic tractability. Nevertheless, some intuition on the spectrum and magnitude of the ionizing 
background during Hell rcionization can be gained by considering idealized cases. We consider two such cases: 1) a 
single quasar at the center of an isolated ionized bubble and 2) a point in a large Hell patch that has yet to be reionized. 

Key insight into the ionizing background is gained by considering relevant physical scales. The left panel of 
Figure |9j compares the mean free paths (calculated as in Appendix [d]) of 1 Ryd (HI ionizing) and 4 Ryd (Hell 
ionizing) photons to the mean separation between the sources of the ionizing background, while the right panel shows 
scales relevant to understanding the possible effects of Hell rcionization by quasars at z ~ 3 — 4. The Hell ionizing 
mean free paths are calculated by converting the HI column densities to Hell assuming a constant Thi = 0.5 x 10~ 12 
s _1 and varying ThcH- Since emission from star- forming galaxies provides most of the hydrogen photoionization rate 
at z > 3 (95J) and the mean separation between L* galaxies is much smaller than the HI ionizing mean free path at 
all redshifts z < 6 considered, 5 it is a good approximation to treat the stellar emissivity as a uniform volume average 
as in equation El It is also similarly the case for HI ionizing quasar emissivity at redshifts z < 3, where quasars are 
relatively abundant and the HI ionizing mean free path large, though with larger fluctuations expected from the 
smaller number of quasars within each mean free path (for more detailed studies of UV background fluctuations, see 
'Zuol|1992bla| |Fardal fc ShiiTT1[T993] |Croft et al.|[TM9l [20021 |Gnedin fc Hamilton|[2002l |Meiksin fc W4utel[2003l |2S04| 
Croft||2004p . Thus, it is a reasonable approximation at all redshifts to calculate the ionizing background between 1 
and 4 Ryd using a volume average emissivity. 



The situation is however quite different beyond 4 Ryd, where continuum opacity owing to Hell dominates. In 
fact, the mean free path of Hell ionizing photons at these energies, which depends on the local Hell photoionization 
rate, is generally smaller than the mean free path of 1 Ryd HI ionizing photons since even quasars produce 

5 At redshifts z > 4, the estimated mean free path relies on an extrapolation of the measured column density distribution and so the 
concl usion should acc ordingly be treated with caution. In particular, the conclusion is likely to break down if HI reionization ends at z ~ 6 



(e.g., |Fan et"alp 002 2006b). 
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Fig. 9. — Important physical scales for cosmological radiative transfer. The left panel compares the mean free paths of 1 Ryd (HI ionizing; 
thick dashed) and 4 Ryd (Hell ionizing; thin dashed) photons to the mean separation between the sources of the ionizing background. 
The different Hell ionizing mean free path curves correspond to different rjj e ii assumed in the calculation and ignore Hell reionization. 



From top down, r Ho ll = W~ 1J , 10~ , 10~ 15 , 10~ lb , and 10~ 17 s _1 , assuming a constant T H i 
between L* galaxies versus redshift is shown by the black points and calculated from measured 
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u ,laxy UV luminosity functions |Stcidcl 

|et al.|1999l|Sawicki fe Thompson|2006||Yoshida et al |2006| |Bouwens et al.|2007||Reddy et al.|2008| lT For the mean spacing betwe en quasars, 
the dotted curves correspond to different lower B- band luminosity cuts and are calculated using the |Hopkins et al.| l |2007[ l luminosity 
function. From bottom up, Lg > 10 40 , 10 42 , 10 44 , 10 46 , and 10 48 erg s . The right panel shows scales relevant to understanding the 
possible effects of Hell reionization by quasars at z ~ 3 — 4. The comoving Hubble radius c(l + z) /H(z) is indicated by the thick solid curve 
and the comoving distance between z = 3 and z = 4 spatial surfaces, labeled the "thickness of Hell reionization" , is shown by the thinner 
solid curve. The dashed curves show the mean free path of high-energy Hell ionizing photons versus redshift assuming a constant ionized 
fraction Xi = 0.5 and that the Hell is uniform distributed, Ai m f p (z) = [o-HcllM^HellC-z)] -1 - The curves, from bottom up, correspond to 
individual frequencies ^ max = (5, 6, 15)vH e ii (sec |7,3.T| for the significance of i^ ma x)- The gray shaded area indicate typical Helll 

ionized bubble radii during Hell reionization. 



relatively few photons above 4 
S = r H i/r H cii = 140 (|Bolton et al.||2006l. 



For example, near the peak of the quasar luminosity function at z = 2.1, 
The relative rarity of quasars and shortness of the Hell ionizing mean 



free path combine to create a situation in which often a single bright quasar contributes to the local Hell ionizing 
flux. This is the case even after Hell reionization has completed and results in substantial fluctuations in the 
ionizing background above 4 Ryd. These fluctuations could be important for metal absorption l ine studies and will 
be addressed in future work (for recent studies of the Hell ionizing background fluctuations, see Bolton et al. |2006| 
Furlanetto|2008 2009 1. Prior to the complete reionization of Hell, ionized bubble walls will further limit the exposure 
of a given point to the Hell ionizing fields of distant quasars. The radii of HcIII bubbles depend, at least until they 
percolate, on the ionizing luminosity of the central quasars and the duration for which these have been shining. This 
gives rise to a wide range of scales, depending on the specific quasar model, but by the middle of Hell reionization 
(determined by an ionized fraction Xi ~ 0.5) bubble radiu s Rb ~ 10 — 100 comoving Mpc (corresponding to a few to 



20-25 proper Mpc at z = 3 — 4) are representative (e.g., Furlanetto & Oh 2008 McQuinn et al. 
average uniform emissivity is then clearly inappropriate. 



2008). A volume 



7.1.1. Quasar Within an Isolated Ionized Bubble 

Consider first a point within an isolated Helll bubble occupied by a single quasar at the center, r = 0. Locally 
neglecting cosmological effects, the specific intensity of a radial ray is given by 



I v = I v ( r = 0)e- T ^ r \ 



(12) 



In addition to the intensity being exponentially suppressed, the spectral shape is altered by the frequency dependence 
of the optical depth: 



J p — T v 

V — = = e ^HcIIe _T "HcIl[' 7 H.Il(f)/'THeIl(^eIl)] „ (g-^Hell ) 



CHell (f)/o , HeII (fcHell ) 



(13) 



where the last equality holds approximately just above the Hell ionization edge and we have neglected the fractionally 
small HI continuum opacity. It follows that the magnitude of the specific intensity is set by the optical depth at the 
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Fig. 10. — Left: Normalized spectral hard ening above the H ell ionization edge. The unhardened quasar spectrum (dotted curve) is 
assumed to be a power law / _1 - 6 (e.g,, |Telfer et al.||2002| . The solid z ~ 2 post Hell reionization curve shows the limit in which 

the mean free path is sufficiently large to contain several quasar sources, with the HI column density distribution and softness parameter 
S = Thi/ThcH = 140 measured at z a 2 ( [Bolton et al.|2006| . This limit is representative of the hardening in the calculations of ij5] The 
dashed curves show the hardened spectrum at different optical depths from a central quasar at the center of an isolated HelH bubble. 
From the thinnest to the thickest curve, t„ HoII = 0.1, 0.5, 1, 2, and 100. The curves are pictorially labeled assuming a Hell ionizing mean 
free path comparable to the bubble size, from close to the quasar to toward the edge of a large bubble and inside a Hell patch. Right: 
Illustration of the absolute effect of Hell attenuation on the full background spectrum. In addition to the spectral harde ning, the spectrum 
is exponentially suppressed above the Hell ionization edge. Recombination reemission is omitted here but discussed in §7.2| 

Hell ionization edge to the source, with the spectral shape entirely determined by the frequency dependence of the 
photoionization cross section at a given optical depth, in addition to the intrinsic spectrum of the source. 



In the left panel of Figure 



10 



we show a quasar spectrum I v {r = 0) oc v is hardened as a function of 
4eII ionization edge from the source. The curves are pictorially labeled assuming a 



tv HoII , the optical depth at the 
Hell ionizing mean free path comparable to the HcIII bubble size, so that tv HoII ~ 1 is near the edge of an isolated 
bubble centered on the quasar. We also show a z ~ 2 post Hell reionization case in which the mean free path is 
sufficiently large to contain several quasar sources, with the H I column density distribution and softness parameter 



S = Tm/Tjieii = 140 measured at z » 2 (Bolto n et al.||2006| ). This limit is representative of the hardening in the 
calculations of f5j Note that as t„ HcII — > oo, the spectral shape can be arbitrarily hardened just above the Hell 
ionization edge, As v — ► oo and cthcIi(^) — ► 0, the spectrum returns to the unfiltered case. 

The rarity of quasars implies that around an individual object the specific intensity obeys equation |12[ in 
which a single source is attenuated with distance, rather than a solution involving a volume average emissivity as 
in equation [2] Why, though, does the ordinary optical depth t v enter in equation [12] instead of the effective optical 
depth f as in equation [2]? Any given light ray is always attenuated according to the intervening ordinary optical depth 
t v . However, the optical depth between two points separated by a fixed distance (at fixed frequency and redshifts) 
fluctuates depending on their particular spatial positions because of the stochastic nature of the intervening absorbers. 
The effective optical depth captures the average attenuation through e~ T = (e _T ). It is an appropriate quantity for 
the ionizing background between 1 and 4 Ryd, where the local intensity is an average over the light received from 
sources in all directions within one mean free path. The radiation above 4 Ryd at a given point in the vicinity of a 
quasar prior to and during Hell reionization will often be dominated by the local quasar and therefore be uniquely 



attenuated as in equation 12 



7.1.2. Point in a Large Hell Patch 

A point within a Hell patch that has not yet been reionized 6 will see a similarly hardened spectrum, but with a 
much stronger suppression at the Hell ionization edge owing to the large intervening optical depth. The optical depth 

6 I n reality, the Hell ionization fronts are quite smooth and extended since they are driven by a hard spectrum (e.g., [McQuinn et al.| 
|2008[ ). Except at the very beginning, few points have been truly untouched by Hell reionization, but the discussion holds wherever the 
ionized fraction has not exceeded, say, ~ 1/2. 
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at the Hell ionization edge, as a function of redshift and path length L, in a medium in which all the helium is assumed 
to be in the form of Hell is given by 

neutral / n ' \ t- r, -, „ I ± ~ * 



CTHcII^HcIlWcIlW-k = 318 



I 



(14) 



4.5 y y 10 comoving Mpc 

In Hell patch, the intensity of the background at the ionization edge is therefore expected be almost entirely suppressed. 
As v — ► oo and CHeii^) - ► 0, however, the optical depth drops quickly and the intensity of the background recovers. 
The corresponding increase of the mean free path with energy leads to the presence of a spatially smooth high-energy 
radiation backg r ound permeating most of the cosmic volume, as seen for example in the numerical simulations of 
McQuinn et al. (20081. The right panel of Figure 10 shows how the fiducial spectrum of the ionizing background at 
z = 3.5, as calculated in ^5| is altered shortward ofthe Hell ionization edge as a function of the intervening optical 
depth. Note, in particular, the tremendous Hell edge suppression even in the moderate case of T UHeI1 — 100. 

7.2. Recombination Lines during Hell Reionization 

The photoionization rate and ionization state of hydrogen are unaffected by the presence of Hell before and during 
Hell reionization apart from a small contribution by photons above 4 Ryd. Consequently, only the Hell recombination 
processes are significantly modified. Of these, the most important is Hell Lya which imprints a distinctive line 
feature at 3 Ryd (Fig. El); Hell LyC and BalC only slightly smooth the spectrum at the Hell and HI ionization edges, 
respectively, and contribute only marginally to the photoionization rates (Figs. [2] and [5]). 



Equations |B7| and |B8| compactly capture the behavior to Hell Lya reemission. As explained in the previous 
section, before Hell reionization begins the background spectrum is almost completely suppressed above 4 Ryd 
by the large optical depth at these energies. Since Hell Lya reemission scales with the Hell ionizing spectrum 
(with saturation in the optically thick limit), it will be absent before the start of Hell reionization. Similarly, no 
Hell Lya should arise within Hell patches during Hell reionization. However, Hell Lya will be reemitted within 
ionized bubbles illuminated by the local quasars. As Hell reionization proceeds, the distance between neighboring 
bubbles should quickly become smaller than the mean free path for the LyC absorption of 3 Ry d ph otons by HI (at 1 
Ryd, Fig. [9] shows the mean free path to be about 200 comoving Mpc at z = 3.5; at 3 Ryd, eq. |D3| predicts the mean 
free path to be longer by a factor of 3 3 ^ -1 ^ « 5) that governs the attenuation of Hell Lya radiation. In the regime 
in which this mean free path contains several ionized bubbles, the volume fraction of reionized Hell can be viewed as 
the fraction of the IGM reemitting in Hell Lya and the Hell Lya reemission line of the background spectrum can be 
expected to be about this fraction times the fully reionized value. 

7.3. Heal Input During Hell Reionization 

The photons that ionize Hell atoms in general carry more energy than the /i^Hcii required. The residual energy is 
converted into kinetic energy of the resulting free electron and HcIII nucleus, with the frequent Coulomb collisions 
leading to rapid thermalization. This process of photoheating is at work at all times and for all species present. Its 
effect is however much more important during reionization, when atoms are being ionized at a much greater rate. 
The effects of Hell reionization on the thermal state of the gas in cosmological simulations has so far generally 
be modeled by artificially boost ing the photoheating rat e calculated from a prescribed spatial ly homogeneous 



background instantaneously(e.g., Bryan & Machacek 2000 Theuns et al. 2002 



Jena et al. 



2005), or ignoring it 



altogether. This approach, of limited physical basis, is a serious limitation of these simulations given the growing body 
of evidence that Hell reionization occurs at observable redshifts z ~ 3—4 and is certain to manifest itself to some extent. 



While cosmological radiative transfer simulations are beginning to self- consistently treat gas thermodynamics 
during Hell reionization (e.g., Paschos et al. 2007 McQuinn et al.| [2008), it is likely that the vast majority of 
simulations performed in the near to moderate future will not explicitly incorporate radiative transfer, either due to 
the computational cost or to the unavailability of an appropriate code. It therefore remains important to develop 
ways of approximately treating the effects of Hell reionization in those simulations. We examine this problem in 
this section. Specifically, we consider the questions: How much does Hell reionization heat the IGM? Over what 
timescale? And how can we approximately model its effects in stand ard cosmological N-body an d hydrody namical 
codes such as GADGET ( |Springel et"al]|2001| |Springel|[2005l ), Hydra flPearce fc Couchman|[l997] , or Enzo ( |0'She a 
et~aI1[2004| ? 



The simple analytic models that follow are motivated b y and owe much to t he physical picture of Hell reion- 
ization suggested by the radiative transfer calculations of McQuinn et al. (20081. We refer to that work for many 
original insights. 

7.3.1. Heat Input Calculation 

In order to gain physical intuition, we begin with a simplified model. Suppose that all the photons up to frequency 
i/ max emitted by a population of sources with intrinsic spectral index ayy are absorbed by Hell atoms. Then the mean 
energy injected into the IGM per ionization is given by 



hv 



Hell 



a uv (1 



x' 



- 1 ) 



auv - 1 (l-£ Q uv) 



1 



Hell 



auv — 1 



(l-auyx^- 1 ), (15) 
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Fig. 11. — Total temperature increase owing to Hell reionization photoheating as a function of the maximum absorbed frequency for 
different spectral indices «tjv °f the reionizing sources. The "expected rang e" corresponds to the bulk of Hell reionization occurring 
between 2 = 3 and z = 4 and is shown only for suggestive purposes. Figure [l2| presents more rigorous calculations that avoid fixing a hard 
cutoff frequency based on the Hell reionization history calculated from the quasar luminosity function. 

where x = ^Hcii/^max and the last equality holds approximation for x <C 1 and auv > 0- This equation neglects 
redshifting of the photons before absorption, which is a reasonable approximation if Hell re ionization lasts Az « 1 

2008 bl |McQuinn et al |[2008|>. Here, the effec ts of 
I'l'ittley fc^eiksin 2007TTBolton et al.|[2008|) are 



at z ~ 3 — 4 (e.g., Furlanetto fc Oh 2008; Faucher-Giguere et al. 



spectral filtering (e.g., Abel & Hachnclt 1999; Bolton et 



al.| |2004 



incorporated in the prescribed frequency cutofi 

If all the helium is initially in the form of Hell and hydrogen is fully ionized, the temperature increase is 
obtained by distributing the injected energy over all particles. After thermal equilibrium has been reached, 



AT, 



Hell 



2 n He 
3fc n to t 



15550 K 



«uv (1 



- 1 ) 



auv - 1 (l-x Q uv) 



1 



31100 K 



/ 0.5 



\auv - 1 



(1 - auv^- 1 ) 



(16) 



Here, n to t = 2nu + 3tih c is the total number density of particles including free electrons. Note that the total number 
of particles is slightly less before Hell reionization owing to the smaller number of free electrons. The fractional 
change of 1/16 is however n egligible. Although the use of a sharp frequency cutoff v max is a simplification of the 
radiative transfer, 
input determined 



McQuinn et al. (2008 1 show that a simple argument like this one gives a good estimate of the heat 
rom detailed radiative transfer simulations. 



What is the relevant value of 



A reasonable guess is the value such that the mean free path of photons 



of this frequency equals the "thickness" of Hell reionization. Photons of higher frequency (and therefore longer mean 
free path) will typically not be absorbed before Hell reionization is complete. The right panel of Figure p\ shows 
where the mean free path intersects the thickness of Hell reionization, assuming that the bulk of the latter takes 
place between z = 3 and z = 4, for different values of ^ max . For this purpose, we calculate the mean free path 
-^mfp(^'max) [nHeiiCHeiU^max)] 1 assuming homogeneously distributed 50% ionized Hell at z = 3.5. Under these 
conditions (accounting for some uncertainty on the thickness of Hell reionization) we expect x" 1 ^8 — 12. Figure 
[TT] shows the corresponding heat input owing to Hell reionization for different value of the spectral index auv- For 
spectral indexes auv ~ 1-5 ([5.1), the heat input depends only weakly on our rough estimate of x^ 1 . 



Having obtained simple estimates for the total heat input during Hell reionization, we proceed to make the 
derivation more rigorous, which also allows us to trace the time evolution of the heat injection. Specifically, we replace 
the sharp frequency cutoff by a calculation taking into account the fraction of photons emitted at each frequency at 
any given redshift that is absorbed during Hell reionization: 

(E )(z) - - ™ - HcIIi v J (17) 



/ z °° dtp dv'/{hv')e^°' com {z')[l - e-^'< *'(*))] 
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where dt — (dz' /c)(dl/dz') and 

t(u', z, z') = £ dz" A„ Ho(z ")[i _ mi ( z ». a uv )]a H cii (u" = ^^^rj) (18) 

is the optical depth encountered by a photon of frequency v 1 emitted at redshift z' before reaching redshift z. Here, 
riHe is the proper number density of helium atoms and a fraction 1 — yni given by the reionization state is assumed 
to be homogeneously distributed in the form of Hell, taken to be the dominant source of opacity. Equation 17 is 
similar to equation |15| but with the high-frequency cutoff replaced by the smoothly varying fraction of photons 
absorbed 1 — e _r for each frequency. In addition, the mean energy injected per ionization is calculated as a function 
of redshift, which allows us to trace the heat input over time. While the homogeneous IGM approximation is 
obviously a simplification, it is a reasonable assumption for this heuristic calculation. In fact, the potential error 
introduced by neglecting the inhomogcncitics is most important for r ~ 1. However, the strong frequency dependence 
of the Hell photoionization cross section implies that the range of photon energy for which r ~ 1 is narrow. 



Morever, for a quasar spectral index auv ~ 1-5, the cruder estimate of Fig. 11 indicates that the heat input is 
only weakly sensitive to the exact maximum energy of the absorbed photons. These effects combine to make the 
uniform IGM approximation relatively robust for this particular calculation. Ultimately, th ough, the calculation i s mo- 
tivated by the fact that it reproduces the results of the full radiative transfer calculations of |McQuinn et al7| ( |2008[ ) well. 

Since at a given redshift z, only a fraction ym of the Hell has been reionized, the temperature increase con- 
tributed by Hell reionization at that redshift, neglecting cooling, is given by 

AT H ciiW = ;|— ym(z)(E t )(z). (19) 
3k n to t 

The temperature of a cosmic gas parcel is in general determined by all the processes by which it gains heat and cools 
as it evolves, including adiabatic heating and cooling, shock heating, photoh eating, Compton cooling off microwave 



background photons, and recombination cooling (e.g., Hui & Gnedin 19971. Instructive intuition can however be 
gained from idealized solutions. 

In the limit of early HI reionization (with the reionization of Hel assumed to proceed simultaneously), the 
temperature at mean density Tq reaches a "thermal asymptote" determined by the competition between adiabatic 
cooling and photoheating and whose value depends on the Hell ionization state. For a power-law background 
spectrum J v cx v~~ ab & just above the ionization edges, a good approximation to the thermal asymptote is given by 

0.53 

T a Sy m P(z) = 2M x 1q4 k(0.464 + 0.536ym)(2 + a^)" 1 / 1 - 7 ( tZ± J (20) 



( Hui fc Haiman|2003[ ) . To first order and ignoring inhomogeneities, the effect of Hell reionization is to inject additional 



heat to each gas parcel. As the universe expands, the extra heat is diluted by adiabatic cooling, T(z) — T(z')[(l 
z)/(\ + z 1 )] 2 , so that an estimate of the overall temperature evolution in the early HI reionization limit accounting for 
Hell reionization heat input is given by 

T(z)^Tr ymp (z)+ [* dz' ^™^ (^j) 2 . (21) 



Figure [12| shows thermal histories calculated using this equation assuming a background spectral index cth s = and 



the Hopkins et al. (20071 quasar luminosity function in the B— band for different spectral indices of the Hell ionizing 
sources auv- Harder spectral indices are seen to result in greater heat injections, which simply owes to the larger 
fraction of ionizations caused by high-energy photons. Moreover, the magnitude of the total heat input as a function 
of spectral index is consistent with the simpler estimates using a sharp frequency cutoff shown in Figure |11| At 
fixed B— band luminosity, harder spectral indices result in higher ionizing photon output rates and thus earlier Hell 
reionization. 



The Helll fraction ym in the above equations is o btained by counting the number of Hell ionizing photons 



emitted by quasars as in Faucher-Giguere et al. ( 2008b I and we have assumed a gas clumping factor C — 5 



7.3.2. Scatter in the Temperature- Density Relation 

The thermal history calculations of the previous section implicitly assumed that the universe is homogeneous at 
a mean density and that Hell reionization happens simultaneously throughout. In reality, the IGM is characterized 
by density fluctuations and the quasars that putatively drive Hell reionization turn on at different times at different 
locations owing to cosmic variance. These inhomogeneities imply that the IGM temperature is not fully described by 
a single redshift dependent number T(z) but in reality exhibits a temperature-density relation T(z; A) with some 
scatter about the mean at each redshift. 
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Fig. 12. — Evolution of the IGM temperature for di ffer ent quasar spectrum indices «tjv in the early HI reionization limit in which the 
temperature would follow the thermal asymptote (eq. p0| with »t,g = in the absence of Hell reionization. Hell reionization is taken to 
occur through the action of quasars, with the B— band luminosity function of |Hopkins et ah] ( |2007| . In the left panel, the evolution of 
the HcIII ionized fraction yui( z ) is artificially fixed to the value calculated for a quasar spectral index «tjv = 1-5- in the right panel, the 
ionization history is calculated consistent with the spectral index of the sources, with harder spectra resulting in higher ionizing photon 
output rates and thus earlier reionization. 



In 



the 

: 2- 



absence of Hell reionization, Hui & Gnedin (19971 showed that the temperature-density relation at 
4 is well approximated by a power law T(z; A) = TqA@ . In the limit of early HI reionization, (3 



z = 'l — 4 is well approximated by a power law l\z; A J = TqA 13 . In the limit of early HI reionization, (3 — > 0.62 as a 
result of the competition between photohcating and adiabatic cooling. We are interested in how this result is modified 
by Hell reionization. Equation [21] can be generalized to 



T(z) 



T a SS ymp (z)A/3 



dz' 



riATHcii(^) f l + z 
dz' V 1 + z' 



(22) 



where (3 is set to the value that would be obtained without Hell reionization and n is a stochastic factor that 
accounts for the fact that different regions are heated at different times by Hell reionization photoheating. Our task 
is then reduced to determining the distribution function of k to estimate the scatter in the temperature-density relation. 



One of the results highlighted by the radiative transfer simulations of McQuinn et al. (12008) is that much of 
the heating during Hell reionization by quasars results from ionizations by the diffuse background of high-energy 
photons with large mean free path s that penetrate int o Hell patches before these are actually reionized by softer 



photons (for a different picture, see Bolton et al.||2008[ ). In this picture, the longer a given region is exposed to the 
high-energy background before it is reionized, the mor e heat it receive s ; regi ons that are reionized last tend to be 
hotter. As an ansatz, again motivated by the work of McQuinn et al. (20081, we may thus posit that k oc i e xp,eff, 



where t e xp,eff is an effective exposure time to the high-energy background. Note, though, that this will not be correct 
at the very beginning of Hell reionization before the background has had time to diffuse. We denote by ZHeii the 
redshift at which a given gas parcel is reionized in Hell and set 



^exp.eff (^Hell) 



dty m (z){\ - yiu{z)). 



(23) 



The effective exposure time is thus the age of the universe at reionization of the gas parcel, weighted by the 
time-dependent ionized fraction, and saturating as the latter reaches order unity. The motivation for the weighting is 
that the heat injection is not only proportional to the raw exposure time, but also to the intensity of the high-energy 
background. The ionized fraction ym counts the number of ionizing photons emitted and is therefore a tracer of this 
high-energy background. The 1 — ym saturation factor approximates the fact the rate of heat input also scales with 
^Heii and is thus suppressed toward the end of reionization. 



The PDF of reionization redshifts is also straightforwardly approximated from the ionized fraction evolution 
since the probability of reionization during a redshift interval scales as the rate at which ionizations occur at that 
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time: 



Energy conservation requires 



-P(zHeii; z) 



yill(ZHell) 1 T7 1 ( 2; HoIl) Z-Hdi > Z 
2 H cII < Z 



(24) 



= 1, so we set k = i xp,eff/(iexp,eff) and in this model the scatter in the temperature- 
density relation ultimately is calculable from the quasar luminosity function. We wish to emphasize that the effective 
exposure time ansatz in equation [23] was obtained heuristically and certainly does not capture the f ull complexity of 
Hell reionization, though it does agree reasonably well with the results of the McQuinn et al. ( 2008[ ) simulations and 
provides a simple way to understand them. 



Figure 13 shows the evolution of the temperature-density relation in this model for the case of a quasar spec- 
tral index auv — 1-5, with = 0. The mean temperature at mean density, (T(z; A — 1)), then traces the 
corresponding curve in Figure U2\ by construction. The slope and scatter of the temperature-density relation are 
however manifest. The underlying slope is set by the early HI reionization limit T — Tq 62 , but the Hell reionization 
heat input somewhat flattens the mean slope. The flattening arises because at low densities that are optically thin 
to the high-energy photons the heat deposition is density-independent in the sense that gas parcels of different 
densities receive the same temperature increment. Since the lower-density el ements are initially co oler, the fractional 
temperature increase is larger for these. In agreement with the simulations of McQuinn et aL] (2008 ) , Hell reionization 
does not produce an isothermal (/? = 0) temperature-density relation. This is, similarly, simply because the initial 
heat in gas parcels above mean density (A > 1) is significant compared to the Hell reionization heat injection and so 
the imprint of this initial heat is not erased. Finally, it is worth recmphasizing that this model predicts a large scatter 
in the temperature-density relation, which may have important consequence for interpreting Lya forest data. In this 
model, though, the scatter arises from the scatter in the Hell reionization times for different gas parcels. Locally, 
neighboring points will be reionized at similar times and we therefore expect the scatter in the temperature-density 
relation to be significantly reduced. 
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Fig. 13. — Evolution of the temperature-density relation in the model of §7.3.2| for the case of a quasar spectral index otjv = 1-5, with 
Q bg = 0- The mean temperature at mean density traces the corresponding curve in Figure [l2| by construction. The slope and scatter of 
the temperature-density relation are however manifest. The underlying slope is set by the early HI reionization limit, T = To A 062 , but 
the Hell reionization heat input somewhat flattens the mean slope. The flattening arises because at low densities that are optically thin 
to the high-energy photons the heat deposition is density-independent in the sense that gas parcels of different densities are heated by the 
same amount. Since the lower-density elements are initially cooler, the fractional temperature increase is larger for these. Hell reionization 
does not produce an isothermal temperature-density relation in this model because the initial heat in gas parcels above mean density is 
significant compared to the Hell reionization heat injection and so the imprint of this initial heat is not erased. 
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TABLE 2 

Photoionization and Photoheating Rates for our Fiducial Model 



z 


r H i a 


T H ei a 


r HeII a 


9Hl a 




*?HcIl a 




AT He II b 




10~ 12 s~l 


10~ 12 s^l 


10~ 12 s~l 


10 -12 eV s"l 


10~ 12 cV s~l 


10~ 12 eV s~ 


1 


T (K) 


0.0 


0.0384 


0.0213 


1.231 x 10~ 4 


0.158 


0.141 


0.0032 




14269 


0.25 


0.0728 


0.0443 


2.956 x 10~ 4 


0.311 


0.299 


0.0073 




14269 


0.5 


0.1295 


0.0860 


6.845 x 10~ 4 


0.569 


0.600 


0.0156 




14269 


0.75 


0.2082 


0.1471 


1.361 x 10~ 3 


0.929 


1.064 


0.0292 
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1.0 


0.3048 


0.2241 


2.317 x 10~ 3 


1.371 


1.677 


0.0476 
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1.25 


0.4074 


0.3076 


3.401 x 10~ 3 


1.841 


2.371 


0.0676 
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0.4975 


0.3843 


4.288 x 10~ 3 


2.260 


3.037 


0.0837 
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0.5630 


0.4446 


4.744 x 10" 3 
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3.579 


0.0920 
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2.0 


0.6013 


0.4856 


4.811 x 10~ 3 
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3.974 


0.0926 




14269 


2.25 


0.6142 


0.5076 


4.511 x 10~ 3 


2.852 
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0.0867 
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2.5 


0.6053 


0.5132 


3.939 x 10~ 3 


2.839 


4.274 


0.0757 
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2.75 


0.5823 


0.5074 


3.223 x 10~ 3 


2.762 


4.246 


0.0622 
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3.0 


0.5503 


0.4942 


2.479 x 10" 3 


2.642 


4.161 


0.0480 




14269 


3.25 


0.5168 


0.4781 


1.812 x 10~ 3 


2.511 


4.047 


0.0351 




11392 


3.5 


0.4849 


0.4617 


1.245 x 10~ 3 


2.384 


3.933 


0.0243 




8007 


3.75 


0.4560 


0.4469 


7.907 x 10~ 4 


2.272 


3.830 


0.0159 




5519 


4.0 


0.4320 


0.4329 


4.818 x 10~ 4 


2.171 


3.737 


9.862 x 10~ 


-3 


3802 


4.25 


0.4105 


0.4203 


2.618 x 10" 4 


2.083 


3.657 


5.639 x 10- 


3 


2619 


4.5 


0.3917 


0.4080 


1.351 x 10~ 4 


2.002 


3.580 


3.133 x 10- 


3 


1770 


4.75 


0.3743 


0.3948 


7.271 x 10~ s 


1.921 


3.493 


1.624 x 10- 


3 


1170 


5.0 


0.3555 


0.3800 


3.724 x 10~ 6 


1.833 


3.393 


8.087 x 10 


-4 


713 


5.25 


0.3362 


0.3647 


2.060 x 10~ 6 


1.745 


3.284 


3.654 x 10 


-4 


394 


5.5 


0.3169 


0.3494 


8.924 x lO" 6 


1.661 


3.169 


1.448 x 10 


-4 


210 


5.75 


0.3001 


0.3327 


5.992 x 10~ 6 


1.573 


3.042 


6.795 x 10 


-5 


112 


6.0 


0.2824 


0.3160 


4.033 x 10~ 6 


1.487 


2.906 


3.461 x 10 


-5 


58 


6.25 


0.2633 


0.2987 


1.657 x 10~ 6 


1.399 


2.762 


1.381 x 10 


-5 


30 


6.5 


0.2447 


0.2801 


9.978 x 10~ 7 


1.305 


2.606 


6.409 x 10 


-6 


15 


6.75 


0.2271 


0.2620 


5.898 x 10~ 7 


1.216 


2.450 


3.022 x 10 


-6 


8 


7.0 


0.2099 


0.2439 


3.430 x 10~ 7 


1.127 


2.292 


1.418 x 10 


-6 


1 



a Optically thin rates. b Cummulative temperature increase owing to Hell reioniz ation. which should be 
incrementally added to gas elements and let cooling to model the heat input (j ]7.3.3| l. 

7.3.3. Implementation in Hydrodynamical Codes 

Cosmological hydrodynamical simulations usually incorporate the effects of a prescribed UV background on 
the thermal histor y of the gas under the assumption of an optically thin plasma (for the relevant equations, see 
Katz et al. 1996a). This assumption manifestly breaks down during Hell reionization and has led simulators to 
artificially increase photoheating rates as a rough approximation of the radiative transfer effects. A more physically 
motivated approach is to increase the temperature of each gas element by an amount (dATn c u(z)/dz)Az (which 
is subsequently let cooling) at each time step Az in the simulation, where AThcIi(z) is pre-computed given the 
desired Hell reionization history as in equation 19 This approach misses the scatter and inhomogeneity of the 



temperature-density relation discussed in the previous section but has the advantage of only requiring an additional 
term in the temperature equation and adding negligible computational overhead while capturing the time scale and 
magnitude of the heat input more realistically. It could conceivably be extended to account for spatial inhomogeneities 
and t he scatter in reioni zation times similarly to semi-analytic models being applied to efficiently model HI reionization 
(e.g., Zahn et al. 20071. An alternative approach would be to replace the optically thin photoionization and heat- 



ing rates by "effective" values that can be substituted into the usual optically thin equations to yield the desired result. 

Table [2] tabulates both the optically thin photoionization and photoheating rates, and the ATneii values to 
use to model the effects of Hell reionization under our prescription for the quasar model with o;qso = 1-6 employed 
in the fiducial background spectrum calculations of ^5] In general, the photoheating rate for species i is given by 



1 

4-7T 



dv 



J v o l {y)(hv - hvi). 



(25) 



The cor responding ionized fraction versus redshift is shown by the C = 5 curve in Fig. 7 of Faucher-Giguere et al. 
(2008bl. In this model, Hell reionization completes by z — 3, with ~ 80% of the ionization occurring between z = 3 
and z = 4. 



COMPARISON WITH PREVIOUS WORK 



As one of the motivations for our calculation of th e ionizing background spectrum was to provide an alternative 
to the widely-used models of Haardt & Madau (1996) (see also Haardt & Madau 20011 and their informally-released 
derivatives, it is useful to directly compare our results with these authors. 
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Fig. 14. — Comparison of our fiducial ionizing background mod el with models inform ally released by Haardt & Madau in 2005 (H&M05). 
Both H&M05 models include a quasar contribution based cm the|Croom et al.| l |2004| luminosity function. One model in addition includes 
a stellar contribution calculated using a |Bruzual &: Charlot| ( |2003p population synthesis code assuming a Salpeter initial mass function, age 
0.5 Gyr, constant star formation, and an escape traction of 107o. First three panels: Detailed spectra at redshifts z = 0, 3, and 6. For this 
comparison, the H&M05 models were normalized by a factor of 0.5 to better match our fiducial calculation. Lower right: Photoionization 
rates of HI and Hell versus redshift. Here, the H&M05 models were not normalized and so the different amplitudes reflect the difference 
in the models as provided. 



Before we do so, we wish to emphasize that in this work we attempted to improve on technical aspects of the 
calculation. In particular, we performed more self-consistent calculations of the ionization structure of individual 
absorbers (fj3| and our treatment of recombination emission (Q was based on approximating the results of these 
photoionization calculations rather than on an escape probability formalism. Moreover, all the numerical calculations 
presented here were performed using an independently-developed code and our empirical constraints (§5| were also 
obtained independently in our previous work. The comparison of the final results against those of Haardt & Madau 
thus provides an indication of the uncertainty in the resulting spectrum. 



In Figure 14 we compare our fiducial model (ignoring Hell reionization) with two models informally released 
by Haardt &c Madau in 20 05 (H&M05; F . Haar dt 2005, private communication). Both models include a quasar 



contribution based on the Croom et al. (2004[) luminosity function. One model in addition includes a stellar 



contribution calculated using a Bruzual & Chariot (120031) population synthesis code assuming a Salpeter initial mass 
function, age 0.5 Gyr, constant star formation, and an escape fraction of 10%. We compare both the detailed spectra 
at redshifts z = 0, 3, and 6 and the integrated photoionization rates of HI and Hell versus redshift. For the detailed 
spectra, the Haardt & Madau models were normalized by a factor of 0.5 to better match our fiducial calculations in 
the ionizi ng range. As a non-negligible uncertainty remains in the amplitude of the intergalactic HI photoionization 
rate (e.g., Bolton et al. 2005 Faucher-Giguere et al. 2008b I, it is fair to renormalize the models before comparing 
them, a procedure which is also fr eque ntly used by simulators to match the observed Lya forest mean transmission 
(see the discussion at the end of ^5.2 1. The photoionization rates shown in the last panel have however not been 



renormalized and so reflect the models as provided 

It is interesting that in spite of the differences in the technical treatment and the independently obtained em- 
pirical constraints, our calculations of the spectral shape in the ionizing range agree quite well with the H&M05 
models between z — and z = 3, suggesting that the calculations are relatively robust in this redshift range. 
The spectra however diverge increasingly toward higher redshifts as a result of the different source prescriptions 
this regime. In our model, the quasar contribution drops more rapidly as z — > oo, while the compensating stellar 
emissivity increases to maintain the nearly flat total hydrogen photoionization rate measured from the Lya forest 
(Faucher-Giguere et al. |2008aTb| . While the |Hopkins et aL fl2007|) lum inosity function we use combines different data 
sets to constrain the faint-end slope up to z = 4.5, thelCroom et al.l (|2004l) luminosity function used in the H&M05 
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models is based solely on the 2QZ survey and is only measured up to z = 2.1. Our calculations are therefore more 
reliable at z > 3. Since we prescribe the HI column density distribution, but self-consistently calculate the Hell 
distribution from the hardness of the background at each redshift, the decreasing Hell to HI ionizing emissivity ratio 
results in a reduction of the Hell to HI ionizing photon mean free path, and therefore amplifies the rHoii/rm evolution. 
This explains the increasingly s tron g Hell break in our model. Although Hell rcionization would modify the results 
above the Hell ionization edge (£ 7.1 1, a prediction of our model is a larger ratio of HI to Hell ionizing flux beyond z « 3. 



A significant difference between our calculations and the original work of Haardt & Madau (19961 is with the 



fraction by w hich recombination emission boosts the photoionization rates. For their original model, |Haardt 
Madau| ( |l996| found rgj th rec /rg£ as high as 1.5 and rg 1 ^ rcc /rg e c n peaking at 1.7. For our fiducial model, we 
found recombinations to be less important for the photoionization rates, with r wlth rec /r rcc < 1.1 at z < 4 for both 
HI and Hell (Figure [5j the boost factor increases somewhat toward higher redshifts as the mean free paths and 
leakage due to redshiftmg below the ionization edges decrease, as discussed in [6.2 1. The analytic model developed in 
Appendix [C] helps us understand our numerical results and giv e us confidence in their accuracy. The differences with 
respect to the boost factors found by iHaardt & Madau ( 1996 1 must partly be due to the different parameters of the 
calculations (e.g., the column density distribution and so urce prescriptions), but m ay also originate from the different 
techniques used. Preliminary investigation suggests that Haardt & Madau ( 1996 ) may have incorrectly used a case B 
recombination coefficient instead of the case A coefficient in a step of their escape probability calculation, resulting 
in an overestimate of the recombination boost which likely explain s a large part of th e discrepancy (F. Haardt 2009, 
private communication). Our results agree better with those of Fardal et al. ( 1998[ ), who find that recombinations 
boost the photoionization rates of both HI and Hell by ~ 20%. 



Recently, Madau & Haardt (2008) proposed a new effect on the spectrum of the ionizing background which 
could be important particularly before Hell is reionized. When a large amount of Hell is present, the opacity arising 
from Hell Ly/3 and higher Lyman-series resonances produces a sawtooth absorption pattern between 3.56 and 4 Ryd 
and the Hell Lya reemission line at 3 Ryd is boosted by resulting degraded photons. We do not include this effect in 
the present work but plan to do so in the future. 

9. DISCUSSION AND CONCLUSION 

In this work, we have revisited the calculation of the UV background spectrum. The three main improvements over 
previous work are: 

• The implementation of new empirical constraints on the sources of radiation based on a d etailed study of 
intergalactic absorption and updated luminosity functions (Faucher-Giguere et al. 2008d|ajb l. In our favored 
fiducial model, star-forming galaxies play a crucial role and a dominate the HI photoionization rate at z > 3. 

• A reexamination of the radiative transfer within individual absorbers and an exploration of the physical depen- 
dences of the calculated background. In particular, we perform more self-consistent photoionization calculations 
including recombination emission, present a new treatment of recombination emission based on them, and clarify 
how the net enhancement of the photoionization rates is influenced by rcdshifting of the recombination photons 
below the ionization edges and their energy distribution. 

• A treatment of the effects of Hell rcionization on background spectrum and on the thermal history of the 
intergalactic medium. 



The main argument supporting a UV background dominated by stellar emission beyond z 



3 is that while the total 
t z = 2 and z = 4.2 (e.g., 



Faucher-Giguere et al. 


2008a|b) 


, the quasar luminosity functior 
•eak and Lya line surveys (e.g., 


i is strongly peaked 


near z=2 (e.g., Hopkins et al.|2007|). 


Thanks to large-scale J 


jyman bi 


jSteidel et al.|1999 


Sawicki & Thompson 2006; Yoshida 


et al.|2006 


Bouwens et al. 2007 


Reddy et al.|2008 


Ouchi et al. 


2008 ) , star-forming 


galaxies are now known to exist nu- 



quasar luminosity function and measurements of Hell-to-HI column density ratios however indicate that quasars do 
contributed a large fraction of a HI photoionization rate are their z ss 2 peak; in our fiducial model, this fraction is 2/3. 

The evidence in favor of a large (and dominant at the highest redshifts) contribution of star-forming galaxies 
to the ionizing background is s upported by related and independent studies. Previous studie s of the HI photoioniza- 
tion rate from the Lya forest (Rauch et al. 1997 Haehnelt et al. |2001| jBolton et al. 20051 have in fact supported 
this conclusion, though with somewhat more leeway owing to larger statistical error bars. The case has also been 
made indepen dently by combining direct measurements of the UV luminosity function of galaxies and of thei r escape 



fraction (e.g., Steidel et al. 



2001 
line 



McDonald & Miralda-Escude 2001 Shapley et al. 



2006; Cowie et al. 2008). Metal 



line studies provide a further line of evidence, indicating that a mix of stars and quasars best fits measured the 
abundance ratios of various ions including CIV, SIV, and OVI ( Boksenberg et aL]|2003| |Schaye et al.|[2003 
et al. 2004 20071. In fact, the ionizing spectrum presented herein could be directly confronted against an 



Aguirre 
further 

constrained by such observations. Theoretical arguments also suggest that star-forming galaxies should dominate 
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early on (e.g., |Springel fc Hernquist||2003| [Hernquist fc Springel|[2003| |Loeb||2008 1 



Although the fiducial model detailed in S|5] fits our observational constraints, it is not at present uniquely de- 
termined. As we explore in fj6] the background depends on the details of both the sources and sinks of radiation. For 
instance, we have adopted a hard spectral index a* = 1 for star-formi ng galaxies between 1 and 4 Ryd based on the 
comparison of theoretical models with observational line diagnostics by |Kewley et al.|(2l)l)l ) . However, in spite of a few 
weak detections, the emergent spectral shape of high-redshi ft galaxies at these energies has yet to b e directly measured 
owing to the large attenuation by the intervening IGM (e.g., Steidel et al.|2001 Shapley et al.|2006 |. Mor eover, present 
popu l ation synthesis models a re at odds with one another in their predictions in this energy range (e.g., Kewley et al. 
|200 1 1 |Boksenberg et aiT||2003[ ), making it difficult to rely on them with confidence. Perhaps the largest uncertainty 
with respect to the absorbers is the abundance and redshift evolution of Lyman limit systems, which determine the 
mean free path of ionizing photons. These are still poorly constrained, especially above z « 4, and introduce a com- 
mensurate uncertainty in the calculation of the amplitude of the ionizing background at these redshifts (e.g., Figure [7]). 

The total ionizing background receives a contribution from recombinations that reemit other ionizing photons. 
We reexamined the physics of this recombination contribution. Interestingly, we find that the enhancement of the 
photoionization rates from recombinations is not simply a function of the number of recombination photons. This 
arises because many of the recombination photons rapidly redshift and leak below their corresponding ionization edges 
as well as from their distribution in energy. Focusing on the HI photoionization rate, the main factors determining 
pwith rccyproc are ^ e parameters of the column density distribution and the LyC recombination line profile. Shorter 
ionizing photon mean free paths, relative to the recombination line width, inhibit leakage. Wider recombination lines, 
associated with higher gas temperatures, on the other hand produce higher-energy photons that receive less weight in 
the photoionization rate owing to the frequency-dependence of the cross section. The steepness of the column density 
distribution also plays a role: for steeper distributions, more of the recombinations occur in optically thin absorbers 
from which essentially all the recombination photons escape into the IGM, in contrast to optically thick absorbers 
that trap a large fraction. In general, r™j th rcc /rgj increases toward higher redshifts as the mean free path diminishes 
and leakage becomes less significant. In Appendix [C] we give a simple analytical model that quantitatively captures 
these effects and agrees well with our full numerical calculations. We conclude that there is not a unique answer for 
how much recombinations boost the ionization rates, but that it depends both on redshift and on the parameters of 
the absorbers and their distribution. 



Our main solutions to the radiative transfer problem assume that the ionizing background is homogeneous. 
This approximation will in particular break down during the rcionization of HI and Hell. As Hell reionization may 



take place under the action of quasars at immediately accessible redshifts z ~ 3 — 4 (e.g., Sokasian et al. 2002 



Wyithe h, Loeb 2003 Furlanetto & Oh 2008 Faucher-Giguere et al. 2008b), we provide a discussion of its effects 
on both the spectrum of the ionizing background and on the thermal history of the IGM. In regions that have yet 
to be reionized, the spectrum is expected to be almost completely suppressed immediately above 54.4 eV by Hell 
absorption. However, the universe remains relatively transparent at higher energies owing the frequency-dependence 
of the photoionization cross section. As the spectrum recovers, a background of > 0.5 keV photons should thus 
permeate the entire universe. 

Another important effect of Hell reionization is to inject heat into the IGM via photohcating. We provide a 
simple analytical model to estimate the overall temperature increase owing to Hell reioni zation based on energ y 
conservation. In this model, which agrees well with the 3-D radiative transfer simulations of McQuinn et al.| ( |2008[) , 
the total temperature increase depends most sensitively on the far UV quasar spectral index. For a value auv = 1-Oj 
the temperature increase could be as much as 15,000 K, though the effect is mitigated by simultaneous adiabatic 
cooling. Harder spectral energy distributions lead to more energy injection. The model is extended to understand 
the effects on the temperature-density relation under the ansatz that the local heat input scales with the effective 
exposure to the high-energy background. The main effects are t o flatten the mean temperature-density relation with 
respect to the early HI reionization limit (Hui & Gnedin 19971 and introduce a large scatter. The flattening arises 



because the high-energy background heat deposition per atom is independent of the local density: the temperature of 
initially cooler lower-density regions is increased by the same amount as that of initially hotter ones of higher density. 
The temperature-density relation does not become fully isothermal as the Hell reionization heat input is comparable 
to the initial thermal energy of the gas and so the memory of the latter is not erased. In our model, the scatter of the 
temperature-density relation originates from the scatter in the reionization times of different gas elements and can be 
calculated from the quasar luminosity function. 

What remains to be done? As outlined above, the ionizing background spectrum is not yet uniquely deter- 
mined and further direct constraints (for example, using metal line ratios) as well as constraints on its sources and 
sinks (which enter the theoretical calculation) are sure to continue to refine it, especially at the highest redshifts. The 
theoretical framework itself needs to be improved to take into account the fluctuations in the ionizing background 
that necessarily e xist at some level and are certain to be i mportant at least during reionization events (for work in 
this direction, see Bolton et al.||2006 Furlanetto 2008 2009]) . In this vein, while our idealized discussion of the effects 
of Hell reionization provides some basic physical understanding, it leaves ample opportunity for improvement. This 
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a particularly exciting area for future progress as studies of Hell reionization are currently blossoming, with much 
HI Lyq fores t data already available and new Hell Lya forest lines of sight that the Cosmi c Origins Spectrograp h 
(Green 20001 to be installed aboard the Hubble Space Telescope is poised to deliver soon (Syphers et al. 2008), 
and the accompanying surge of theoretical interest. In particular, severa l groups are now beginning to nu merically 



Paschos et al. 2007 



tackle the full problem of 3-D radiative transfer and its thermal effects ( pokasian et al.||2002 i 
[McQ uinn et al. 2008; Mascl li et al.|2008| . Another interesting observational opportunity for probing the extragalactic 
radiation background is provided by the newly-launched Fermi Gamma-Ray Space Telescope. Indeed, the intervening 
extragalactic background light should attenuate the 7— rays from distant sources through electron-pos itron pair 
production and thus give us an additional ha ndle on it, particularly below the HI ionization edge (e.g., Madau & 
Phinney|[l996| |Qh]|2001a| |Razzaque et~al~]p008 1 . 



To facilitate use and extension of the results presented in this work and their comparison with observations, 
numerical tables are available in electronic form on the web at |http://www.cfa.harvard.edu/~cgiguere/uvbkg.html| 
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APPENDIX 
A. PHOTOIONIZATION CALCULATIONS 

In this Appendix, we detail the code used to calculate the photoionization equilibrium structure of individual cosmic 
absorbers outlined in ^3] 

We approximate the absorbers as sheets infinite in extent but finite in thickness, with geometry defined in 
Figure [15] Although this geometry is restrictive, the calculation is otherwise three-dimensional in the sense that it 
takes into account that rays incident at different angles encounter different optical depths. This geometry is clearly 
adequate for the sheets of the cosmic web, but somewhat inexact for the filamentary and clumpy structures, and the 
results may therefore be off by a corresponding geometrical factor. Nevertheless, the calculations retain the essence 
of the problem and significantl y improve over previo us work that either as sumed a sem i-infinite geometry and an 
escape probability formalism ( |Haardt fc Madau| [1996) or gray cross sections ( |Fardal et al.||1998| . We also for the first 
time self-consistcntly treat the coupling between hydrogen and helium arising from their recombination emission, as 
explained below. A more accurate approach would consider a distribution of three-dimensional absorber geometries 
(e.g., obtained from a cosmological simulation) for each line-of-sight optical depth considered, but would be much 
more involved and is beyond the scope of this work. 



We assume th at the slab is com posed of hydrogen and helium, with cosmic mass fractions X — 0.75 and 

x 10 4 K, as estimated for optically 
see also Schaye et al. 2000 and Ricotti et 



al 

with iVjji > 



Buries et al. 


2001). The temperature of the gas is set to T = 


absorbers (M 


LcDonald et al. 


2001 


Zaldarriaga et al. 2001 


see 



2000). While the gas temperature may differ and be more complex m structure in dam ped Lya absorbers (DLA) 



2 x 10 cm that are able to cool and form stars (e.g., 



117.2 



cm 
slab 



Wolfe et al. 
1 are comp 
is assumed 



2005J, the detailed properties of 
etely opaque to ionizing photons 
to be equal to its Jeans length 



these absorbers are not crucial since absorbers with 2Vhi ^10 
regardless of their exact column density. The thickness of the 

Lj = yJir^fkT /Gpiim Vl 7 which is both theoretically motivated and provides a good match to observations ( |Schaye| 
2001). In practice, we prescribe the physical mass density p and derive the ionic column densities in phot oioniz at ion 
equilibrium so that we do not need to explicitly assume a relation between Nm and nn- Because pressure also 
smooths the gas on this scale, the absorbers are assumed to be uniform in density, but we have explored other 
density profiles and found that our numerical results are only marginally affected, and our broad conclusions unaltered. 



Specifically, we solve the following set of equilibrium equations under external illumination from both sides by 
an isotropic radiation field of specific intensity equal to the cosmological J v at each point in the absorber: 

a H i(T)n H ii7i = r H in H i, (Al) 

ctHci(T)n HcII n c = r Hc i^Hoi, and (A2) 

aHeii(J 1 )n H oiii"-c = r Hc ii«Hcii (A3) 

7 Here, 7 is the adiabatic index and the mean molecular weight of the gas. These are taken to be 5/3 and 0.59, respectively, 
corresponding to a monatomic and fully ionized gas of cosmic composition. 
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subject to the constraints nu = ^hi + nmi, "-He = ^Hci + «Heii + ^Hciii, and n = nnn + HHeii + 2nHein- The 
photoionization equilibrium assumption is generally accurate as the ionization timescale is much smaller than the 



Hubble time and collisional processes arc negligible at the densities and temperatures considered (Haardt & Madau 



1996) 



To close the system of equations |A1| - |A3[ we must specify how to calculate the photoionization rates. Let 
J v = (1/4tt) J (011^,(8) be the angle average of the total specific intensity at any given point. Then the local 
photoionization rate for species i £ {HI, Hcl, Hell} is given by 

/■OO T 

r i (x)=47r/ ^-ai{y)J v {x). (A4) 

The specific intensity along a particular ray will in general depend on the angle 6 of incidence inside the slab because 
different rays encounter different optical depths. We take the total specific intensity at any given point to be the sum 
of the radiation originating from the external background and of the radiation from the recombination processes to 
HI and Hell (HI LyC, Hell BalC, Hell Lya, and Hell LyC) within the absorber calculated as in jj4]and Appendix [b| 
For this purpose, the recombination line photons are treated as 8— functions, which is a good approximation within 
individual absorbers. Each component is attenuated in magnitude depending on the optical depth from its source 
following the usual transfer equation. 

In particular, our approach explicitly takes into account the coupling between different species arising from 
the reabsorption of recombination photons by a different species. For example, ionizing Hell LyC recombination 
photons will in general not only be reabsorbed by Hell, but also by HI and Hel. Be c ause of this explicit treatment 
of recombinations, case A recombination coefficients are appropriate in equations |A1| — A3 As Hel is found to play a 
negligible role in our calculations, we do not treat its recombination processes in detail 
effect by using the case B recombination coefficient in the optically thick regime. 



but we do approximate their 



Our calculations use 100 spatial bins, 100 logarithmically-spaced frequency bins, and 20 angular bins covering 
< 9 < 7r < 2. The equations are solved iteratively until convergence to better than one part in 10 3 is attained at 
each point. 
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Fig. 15. — Definition of the geometry for our photoionization calculations. Absorbers are modeled as sheets infinite in extent, but finite in 
thickness. The ionization state of the ions of hydrogen and helium are calculated for a prescribed gas temperature T and isotropic external 
radiation field of specific intensity . 



B. RECOMBINATION EMISSION 



Here we provide the technical details of our treatment of recombination emission from individual absorbers outlined 
in f|4j including the analytic approximations to the self-consistent numerical photoionization results. 
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General Formalism 

Consider a generic line recombination process with source function S£ ec = j u ec /a u , where 

jl cc = kl ^a rcc (T)n l+1 n e <fi 1 . cc (v) (Bl) 

is the relevant emission coefficient and a v = 'Yl,i n i' J A 1 A ^ s the absorption coefficient accounting for absorption by 
all species. Here, v rcc is the frequency of the recombination line of interest; we make the approximation that the 
line is sufficiently narrow that it can be represented by a single frequency for ene rget ic purposes and the frequency 
dependence of the line profile is captured by the function (jW^), discussed in S: 4.2 The effective recombination 



coefficient a ICC (T) accounts for all the channels leading to the transition of interest. We use the subscript "i + 1" to 
refer to singly ionized species i. For instance, if i is HI, then i + 1 is HII, and if i is Hell, then i + 1 is HcIII. 

A given ray crossing a slab absorber with incidence angle 9 w ill emerge with extra rec ombination photons 



following the general solution to the radiative transfer equation (e.g., Rybicki fc Lightman||1979 1: 



pT v (X — OO) 

irWn; x = oo, M )= / dir'J^e-^-^^STirl), (B2) 



where t v is the optical depth normal to the slab and /i = cos 9. The recombination intensity can also be calculated in 
this way at any point x within the slab to model the effects of recombinations on the ionization structure of the slab 
itself (93^ and Appendix |A[) . Since d 2 N/dzdNm is the observed column density distribution, the appropriate value of 
/i to use in calculating the cosmological emissivity (eq. [8]) is one and we define Il ec (Ni) = I u cc {Ni, x = oo, /i = 1). 

In the optically thin limit t v (x = oo) <C 1, 

!v (N m ) NiT t (j) ICC (is), (B3) 

so that the recombination emission is proportional to the column density of the absorber (Ni), times the number 
of incident ionizing photons (r^), times the fraction of recombinations that lead to the recombination process of 
interest (a lec (T) / af (T)) . As an absorber becomes optically thick, the number of incident photons it absorbs saturates 
as 1 — e~ T " , so that we expect that the recombination emission will also saturate accordingly. As the numerical 
calculations in Figures [2] and [3] show, this is the case and will be the basis for the analytical approximations to the 
reemission that we develop next. 

Analytic Approximations 

A good analytic approximation to the emergent reemission from an absorber, capturing the optically thin limit and 
the optically thick saturation, is given by 

(T) 

4tt a?(T) 



I^(Ni) = ,1 CC ^J mmiNi, 2V,, thresh )r^ re c(^), (B4) 



where and iV t hrcsh is a threshold column density, near the optically thick transition, at which the recombination 
intensity saturates. To make the approximation smoothly varying with column density, it is convenient to make the 
replacement min(A r ,;, -/V^thrcsh) — * -/Vj !t hresh(l — e^ Ni ^ Ni - thIOBh ), which we do throughout in our numerical evaluations. 
By inspection, we find that TVhi, thresh = 10 16,75 cm~ 2 and iVneii, thresh = 10 17 3 cm -2 give good approximations to 
the full numerical results for the HI LyC arid Hell LyC processes, respectively (Fig. 2). 

For Hell BalC reemission, a more robust approximation is obtained by noting that in the optically thick regime, the 
absorbers have a "skin" (analogous to a region behind an ionizing front) that is nearly uniform in HI and Hell. While 
the recombinations are to Hell, the opacity to the 1 Ryd recombination photons owes to HI. Because this skin is 
optically thick at the recombination energy, the emergent intensity simply approaches the source function. The value 
of the source function near the outer edge of the absorber can be approximated using optically thin photoionization 
abundance ratios: 



rHell BalC / at \ ^ e HeII BalC / e ,i-^ ^ ^^HoII BalC aHcILn=2(r) ^thinrifcll / / \ /T>t;\ 

i v (JVhi) -> b y (Skm) w ^ — — r-<Pr CC \v). (B5) 

For this process, we approximate the reemission as the minimum of the optically thin limit and this optically thick 
result. 

As Lya is a resonant transition, its recombination photons will scatter until they diffuse out of the absorber 
unless they are destroyed by dust, metals, or HI continuum opacity in the case of Hell Lya which can ionize it. 
Of particular interest for the destr uction by metals is the coincidence between the OIII A303.799 line and the Hell 
Lya line at 303.783 A ( Bowen|l934 1. Studies of this process (including HI opacity) in the context of dense and enriched 
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plane tary nebulae and galactic nuclei ( |Weymann fc Williams|[T969 Kallman fc McCray| [l980 Eastman & MacAlpine 

naiority ofThc Hell Lya emission created by recombinations ditfi 



1985 1 suggest that the vast majority of the Hell Lya emission created by recombinations diffuse unimpeded into the 
1GM for Lya forest systems. Only the most metal-rich Lyman limit systems and damped Lya absorbers potentially 
lose a significant fraction of their Hell Lya recombination photons to OIII, which should have little global impact. 
Similarly, only the highest column density and most chemically evolved systems are likely to contain enough dust to 
efficiently destroy HI Lya photons. We therefore adopt a simplified treatment of HI and Hell Lya recombination 
emission in which we assume that all recombination photons escape into the IGM. Instead of integrating the source 

function over the absorber, we thus simply integrate the emission coefficient j™/ HcI1 Lya. 

/oo 
dx jRi/nen Lya^ (B6) 
-oo 

By the preceding arguments, this is likely a reasonable approximation, but certainly an upper bound. We assume 
case B conditions in which higher Lyman-series photons are ultimately degraded into lower-energy photons, with the 
ultimate production of extra Lya photons, which we implement by taking the appropriate case B emission coefficient. 



For the Lya processes, particularly good approximations to the numerical results can be developed as the op- 
tically thick limit can be accurately estimated. The key observation is that the fraction of recombinations ultimately 
leading to the reemission of a Lya photon is a fixed number (at a given temperature) equal to the sum of all the 
production channels allowed by the selection rules. Recombinations directly to n = 1 produce LyC photons, but 
unless the recombination occurs near the skin of the absorber this photon will be reabsorbed before escaping and one 
must account for the probability that it will ultimately result in a Lya photon. Averaging over angles, the number of 
ionizing photons incident from one side of the absorber that are absorbed per unit time per unit area per unit solid 
angle is given by 

7>1 7>OG 7 />00 7 

K hs {Nm) = J o dp I _j-[i_ e WM ] = ^ ±_j ni _ Tur{ _ lt Tv)]) (B7) 

where here T(.,.) is the upper incomplete gamma-function and is unrelated to the photoionization rate, and i £ 
{HI, Hell} indicates the species of interest. Since in equilibrium the total number of recombinations equals the 
number of ionizations, the Lya reemission in the very optically thick regime can then be written as 

^m/Hen Lya (Nm) = ^j^n **°( Nm)S ( v , ffl/HcII LyQ ), (B8) 

where /^ hlck = 0.4, set by matching the numerical calculations (Fig. [3]), accounts for the efficiency of Lya photon 
production and geometrical effects. In general, we again approximate the reemission as the minimum between the 
optically thin and optically thick results. 



Note that an alternative approach for approximating the recombination emission can be derived from the fact 
that the cosmological emissivity (eq. [8]) ultimately does not depend on the detailed function J™ C (A^) but on its 
integral over the column density distribution. The pa rt of the integral inv olving Ni can thus be factored out and 
pre-computed for a given column density distributio n. [Fardal et ah ( 1998 1 used an approach along these lines and 
we also make use of this fact in the next section (eqs C5 and C7 ) in developing an analytic model for the integrated 
recombination contribution to the ionizing background 



C. ANALYTIC MODEL FOR THE RECOMBINATION CONTRIBUTION TO THE IONIZING BACKGROUND 

Here we provide a quantitative analytic basis for und ersta nding for the full numerical results on the contribution of 
recombinations to the ionizing background obtained in £6.2 We focus on the HI photoionization rate and assume that 
the mean free path is sufficiently short that the local source approximation of equation |D1| applies, which is valid at 
z > 2. Based on fjljand in particular Figure [2] we assume that HI LyC reemission is the dominant process. Following 
the notation of the main text, let us define r*j th roc = r no rcc 
rate that is contributed by recombinations. 



HI 



r^f so that Tjjj is the portion of the photoionization 



As outlined in |6.2[ the recombination contribution depends on several factors. Quantitatively, we write 

pre 
1 u 



prcc 

1 HI 
"pwith rcc 
1 HI 



HI 

"pwith rcc 
1 HI 



X /leak X / CTJ 



(CI) 



where the first term is the maximum that would be attained if all the reemitted photons contributed to the ionizing 
background and all had frequency z/jji, /leak accounts for leakage of the photons below the ionization edge, and / CTi 
accounts for the frequency-dependence of the photoionization cross section and recombination line profile. 



We first consider the maximum contribution by assuming that a negligible fraction of the recombination pho- 
tons leak below the ionization edge before being reabsorbed and taking the recombination line profile to be a pure 
(5-function, e™ c = e^ c ,<5(^ — v m )- Then, 

r rcc , f°° dv Tree , s A ^P (>Hl)gHjjg , . 

F H i = 4* / Vy J v a m {v) = j— . (C2) 
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Combining equations [8] and B4 for the cosmological recombination emissivity and the special case of HI LyC gives 

, dz 
dl 



~rcc 

fHI 



&HI,n=l(T) witll rcc ( 

1 T- 



HI 



00 d 2 N 



-JVHl/JV m ,thr 



and therefore 



where 



t^icc 
1 HI 



T^with rcc 
1 HI 



/off' = Al mip {v m )a m — J 



dzdN m 
d 2 N 



eff - 



dNm ^ n , r — ^HI, thresh (1 



-JVhi/JVh 



.thresh^ 



(C3) 
(C4) 
(C5) 



is a dimensionless efficiency factor whose value depends on how recombi natio n emission saturates with column density 
relative to absorption (the mean free path term). Using equations [3] and D2 we can express the mean free path at the 
ionization edge as 

dl r " x 



A^mfp^ffl) 



dz 



d 2 N 
dN m (1 
ozoN m 



-TfflJVffl 



and thus 



/eff 



/ °° dN m d 2 N/dzdN m (a Hl N HItthlcsh )(l - e -^H./Af H i,th„.h) 
J °° dN m d 2 N/dzdN m (l - e -°™ N ™) 



(C6) 



(C7) 



For a power-law column density distribution as in equation [5j the redshift dependence cancels out in the ratio and the 
integrals over iVjji can be done analytically in terms of gamma-functions. The latter also cancel out, leaving a very 
simple result: 



/off — (CHICHI, thresh) 



1-p 



In we found that Nm 



thresh 



10 



16.75 



cm 2 provides a good approximation for HI LyC. Since cthi^hi 



thresh 



(C8) 
< 1, 



equation |C8| shows quantitatively how the r ecom bination contribution increases with the steepness of the column 
density distribution (large f3). As noted in {6.2 this physically arises because a larger fraction of recombinations 
occur in optically thin absorbers which return more of their recombination photons into the IGM. 

For a hydrogenic atom of integer charge Z, the fraction of recombinations that are to level n is 

a H i,n(T) n- 3 e^/ kT Ei^ n /kT) 



where <pi — Z 2 hvm/i 2 is the ionization energy of level i and Ei{x) = dx'e~ x jx' ( Cillie|[l932 ) . 



(C9) 



The actual contribution of recombinations to the ionizing background is smaller because many recombination photons 
rapidly redshift below the ionization threshold. What fraction of photons are lost through this leakage? Focusing 
again on LyC recombinations, a photon is reemitted just above the ionization edge with energy v' = vm + Av, where 
Avjvm <C 1. Supposing this photon is reemitted at redshift z roc , it will redshift below vm after traveling a proper 
distance ^ 

AJ leak (z r c C ; Av) K. — (z rec ) (1 + Z r cc)- (CIO) 

az fjji 



The local source approximation of equation |D1| is valid because in this regime we can write J„ = 
(47r)~ 1 Jq 00 <iZe„e~'/ A ' mtp , with the emissivity treated as a constant. The leakage of recombination photons owing 
to redshifting implies that the emissivity should really integrated over a maximum distance Al\ c& \: 

JT = ^ j"** dK cc e~ l / Al ^ = i- A/ mfp e r ; c [1 - / lost ] , (Cll) 

where /i os t(zrcc; Av) = e - A; i<=ak/A; mfp j g f rac tion of recombination photons emitted with frequency above the 
ionization edge that are lost to redshifting. 

The overall fraction of photons lost is then an average over the recombination line profile: 

/•CO 

/lost(-Srec) = / ^'0rec(^')/lost(z r ec; Av = v' - V m ) . (C12) 



The result is simplified if the line profile is taken to be purely exponential, 4> Tec (v') = (h/kT)e hA "/ kT for Av > 0, in 
which case 

/lost^rcc) = — -; —j — , (C13) 

1 + y{z Tec ; A< mfp , T) 
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Fig. 16. — Dependences of the recombination contribution to HI photoionization rate calculated with the analytic model of Ap- 
pendix [C] Shown are ratios of the total HI photoionization rate, including recombination emission, to the same calculation without 
recombination emission. We vary the parameters of the column density distribution and temperature of the absorbers. Top left: 
C = 0.125, 0.25, 0.5, 0.75, and 1 from bottom up. Top right: (3 = 1.2, 1.3, 1.4, 1.5, 1.6 and 1.7 from bottom up. Bottom left: 
T = 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 X 10 4 K from bottom up. Bottom right: 7 = 1.0, 1.25, 1.5, 1.75, and 2.0 from bottom up. The 
first three panels can be directly compared with the corresponding full numerical calculations shown in Figure J8] The main features and 
dependences are well reproduced by this simple model. The differences at z < 2, where the analytic model is inapplicable, arise because 
the local source approximation used becomes invalid as the radiation field begins to be limited by the cosmological horizon. 

where 

y(z Tec ; AU P , T) = r (C14) 

Al m{p {dz/dl)/{l + z rcc ) 

The mean free path depends on redshift and on the parameters of the column density distribution. Using the analytical 



expression D3 for Al m f p , we can express /i ost directly in terms of these basic parameters: 



/lost^rec! C, 7, f3, T) — ^rp/ . U+ „ T/ , (C15) 

1 + CT(2 - j3) (1 + z rcc ) 1 +TfcT//i^Hi 



and finally define / loak = 1 - /] 



lost ■ 



The last term to consider is the suppression factor that arises because the photoionization cross section enter- 
ing in the photoionization rate from recombinations is frequency-dependent and photons are really reemitted with 
finite energy above the ionization threshold. This is simply a frequency average of the cross section, over its maximum 
at Vm: 



.L - J°° <h/<t> ne {i/)( V / m )- i = e-^/ fcT r(-2, /,7',,,/A-n. (Cl(>) 



where we have again approximated the recombination radiation to have a purely exponential profile and that the 
cross section scales as z/ -3 just above the ionization edge. 



In Figure 16 we combine these analytic results and show how rgj th rec /rg^ rec = 1/(1 - r^f/rgf 1 rec ) varies 
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as a function of the parameters of the column density distribution and the temperature of the absorbers. Note 
that the main features and dependences of the corresponding full numerical calculations shown in Figure [8] are well 
reproduced by this simple model. The differences at z < 2 arise because the local source approximation used in the 
analytic calculations becomes invalid as the radiation field begins to be limited by the cosmological horizon; the 
"effective mean free path" is then shorter than the mean absorption distance and more recombination photons are 
thus in reality retained in the ionizing range. 

D. SPECTRAL FILTERING 

The spectral index of a radiation background in general differs from the spectral index of its sources owing to 
filtering along the line of sight. In this section, we explore different filtering cases relevant to the ionizing background 
in order to provide physical understanding of our numerical solutions of the radiative transfer equation and provide 
analytical results referred to in the main text. 

An important characteristic of the ionizing background at high redshifts is that it is local in the sense that 
the specific intensity J„ depends only on the local value of the specific emissivity e v : 

M*) ~ T^AZ mfp (^ z)e v (z). (Dl) 

47T 

Here, A/ m f p (j/, z) is the mean free path of photons of frequency v and redshift z and is given by 

AU>K = - )"W (° 2 ) 

This limit of equation[2]is valid whenever photons are absorbed so close to their point of emission that they redshift only 
negligibly. For HI ionizing p hotons of waveleng th 912 A, the "breakthrough" point above which this approximation 
holds is approximately z = 2 (Madau et al.|1999|. In this regime, cosmological effects are unimportant and calculations 



can be performed in ordinary Euclidean geometry. Then the effective optical depth at frequency v over a proper length 
I at redshift z can be written as f(y, z, I) and AZ m f p (^, z) — (df/dl) -1 . A useful analytical expression for the mean 
free path can be obtained at frequencies where the photoionization cross section o~i oc and the column density 
distribution is described by single power laws (3 and 7 (see eq. [5]). In the case of HI, 

a/ ( ^~ hV' M) 1 r™ 

AUp( '°' Za) rtM^f UJ (i + ^H( Zo y (D3) 

and it is straightforwardly generalized to any other single ion. 



An implicit assumption in equation Dl (as well as throughout much of this paper, such as as in eqs. [l] and 



Pj) is that a sphere of radius one mean Tree path contains sufficiently many sources that their effect is well-captured by 
the use of a volume-averaged uniform emissivity. This is generally valid when the sources are star-forming galaxies, 
which are very numerous. Quasars, however, are much rarer and this assumption in general fails. (See Figure [9] and 



the discussion in {7.1 ) We will therefore also consider the case of filtering of radiation from an isolated source. For 
each case, we will consider both the cases in which the intervening absorbers are uniformly distributed and the one in 
which discrete absorbers are Poisson-distributed following a column density distribution. 

Uniform Emissivity and Absorbing Material 
In the case in which both the emissivity and the absorbing material are spatially uniform, 

JA Z ) ~ T-Mnfp^, z)e u (z) with f(v, z. I) = n ahs a ahs (iy)l, so that J u (z) « ^- €v ^ Z (D4) 

47T 47T ni<7i{V) 

where rii is the number density of the absorbing material and <Ji(v) is its cross section. For photons with ^hi < v < i^Heii, 
the dominant source of continuum opacity owes to HI photoionization, so rij = rtui and Oi(y) = cthi(^)- For v > i^Heiij 
HI continuum opacity is fractionally small in the cosmological context and we can take ni = nneii and o~i{y) = OKeiiiy). 
Since both absorbing ions are hydrogenic, they similarly harden spectra following J„ cx e v /oi(v) e„f 3 . The last 
equality is approximately valid near the ionization edge. At high energies, the photoionization cross section decreases 
more slowly and the hardening becomes negligible. 

Uniform Emissivity and Discrete Absorbers 

If the emissivity is spatially uniform, but the absorbers are discrete and Poisson-distributed, the result is similar, 
but with the effectively optical depth calculated from the column density distribution as in equation [3] If the column 
density distribution of the absorbers follow a power-law dN/dNi cx N~ ^, then just above the ionization edge where 
the cross section <Ti{v) oc f -3 , the mean free path A/ m f p (/;, z) oc i> J ^ _1) 



2008b I and therefore J v oc e^zr 3 ^ 1 ). 



Zuo & Phinney||l993 Faucher-Gi 


guere et al. 


measured by |Misawa et al. 


(2007 


), in which 



case the hardening a — > a — 1.2 is significantly weaker than the a — > a 
above the HI ionization edge. If the proportionality factor rj between Nm and -/VhoII were a constant throughout, the 
radiation above the Hel l io nization edge would be identically hardened. However, the complex behavior of 77 in the 



optically thick regime 3.2 ) may alter this result; this behavior is taken into account in our numerical calculations. 
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Point Source and Uniform Absorbing Material 

The case of an isolated point source is quite different. In this case, the specific intensity (which can no longer be 
assumed to be isotropic, so we denote it by I v instead of J„), is exponentially suppressed with increasing optical depth 
from the source: 

I v = h{l = 0)e- T "W . (D5) 

We again assume that HI continuum opacity dominates for ^hi < v < i/Heii and that Hell continuum opacity dominates 
for v > fHcii- In each regime, t v (V) — niOi(y)l. The spectrum is attenuated relative to its value at v\ by a factor 

k. = = e ^e-nK(-)/-i(^)] K (g-r^MMM ra ( e -r„ ( )(«Mr 3 _ (D6) 

where the last equality again holds approximately just above the ionization edge. 

Point Source and Discrete Absorbers 

For a point source attenuated by Poisson-distributed discrete absorbers, we simply replace the ordinary optical depth 
by the effective optical depth to obtain the average spectrum: 

/„(/) = = 0) e -^ l l (D7) 

Drawing on our previous results, 

k = = e**e-*" oc e-- « {e^f M ^ . (D8) 



It is worth noting, as explained in §7.1.1[ that at fixed optical depth from a source the hardening is the same 
regardless of how (smoothly or discretely) the intervening material is distributed. The different average hardening 
differs in the discrete case really because of the added stochastic nature of the intervening optical depth at fixed 
distance. 



E. ATOMIC PHYSICS 



The radiative transfer calculations in this paper are ultimately roote d in atomic physics. T he recombination rates 



and Hcl photoionization cross sections are taken from the appendix of |Hui fc Gnedin (19971. For the Lya emission 
coeffic ients and the HI and Hell photoionization cross sections we take the expressions given by Osterbrock &; Ferland 



(20061. In particular, given the prominent role it plays in our calculations, the photoionization cross section of a 
hydrogenic atom of atomic number Z (Z = 1 and 2 for HI and Hell) is given by 

A ^4 exp {4- [4tan2MA} (vU 
Z 2 - \ v I 1 — exp(— 2n/e) 

for v > V\ and otherwise. Here, 

2 9 7T 

A Q = —^airal = 6.30 x 10" 18 cm 2 , (E2) 

where a is the fine structure constant and ao is the Bohr radius, e = dv ~~fv\ — 1, and hv\ = Z 2 hvY{\ = 13.6Z 2 eV. 
Just above the photoionization edge v\, i.e. for e -C 1, <Ji{y) oc is~~ 3 but the cross section drops more slowly as v — > oo, 
explaining the lack of spectral hardening in this limit. 

Also of interest is the recombination line profile for free-bound transitions, which is important in determining 
the contribution of recombination emission to the ionizing background ([Q. The probability that a recombination 
yields a continuum photon of frequency between v and v + dv 7 by definition, is (j)- cec (y)dv . The velocity u of the 
recombining electron relative to the nucleus 8 is related to the frequency by hv — m e u 2 /2 + hv rcc , where v rcc is the 
ionization edge frequency. Now, the probability that the recombining electron has velocity between u and u + du 
scales as the probability that an electron in this velocity range recombines when it encounters a nucleus, times the 
number of electrons with velocity in this range. If <7 rcc (u) is the velocity-dependent recombination cross section, 
then the first term is oc er rec (u)u. In thermal equilibrium, the second term is given by the Maxwell-Boltzmann speed 
distribution, /m-b(w) oc u 2 exp(— m e u 2 /2kT). Thus, 

<t> xec (v)dv oc a Tec (u)nfu-n{u)du. (E3) 

The Milne detailed balance relation relates the recombination cross section to the corresponding photoionization cross 
section: 

a lec (u) = 2 ( ] (Tiiy) (E4) 

\m e uc J 

8 Note that since m p /m e PS 2000, the nuclei can be assumed to be at rest in the frame of the gas. 
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(e.g., Osterbrock fc Ferland|p006 1 . Using the approximation Oi{y) cx v 3 , we can solve and find, after normalizing, 

, , , {v/v^)' 1 exp(-hv/kT) 6(i> - v TCC ) 

4>^c(v) = — — — . E5 

T(0, hvr CC /kT) v Tec 

This profile will also be thermal broadened and shifted owing to the peculiar velocity of the emitting gas; these 
corrections are however negligible in comparison to the width of the recombination line profile. For instance, the 
recombination line width Ais rcc /v rcc ss kT/hv rcc 0.13 for HI LyC at T = 2 x 10 4 K, while the Doppler broadening 
at the same temperature Aisd/v tcc ~ 0.002 and the peculiar velocity shift Ai/ pcc /^ rcc w 0.0003 for v pcc = 100 km s . 
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